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Abstract 


In  this  report,  an  experimental  technique  is  described, 
by  which  some  basic  properties  of  viscoelastic  materials  such  as 
phase  velocity,  loss  factor  and  stress  relaxation  function  could 
be  derived.  The  technique  is  based  on  the  propagation  and  reflec¬ 
tion  of  longitudinal  stress  pulses  in  a  system  of  two  bars,  one 
of  which  is  elastic  with  known  properties,  and  the  other  visco¬ 
elastic  with  properties  to  be  determined.  It  is  shown  that  the 
unknown  properties  of  the  viscoelastic  bar  may  be  determined  from 
the  shapes  of  the  incident  and  reflected  stress  pulses  which 
travel  In  the  elastic  bar.  It  turns  out , however / that  the  only 
material  property  which  may  be  derived  with  reasonable  accuracy 
is  the  phase  velocity  for  the  basic  frequency  of  the  pulses 
propagating  along  the  bars.  All  the  other  material  properties 
proved  to  be  extremely  sensitive  to  small  experimental  errors 
and  hence  no  reliable  values  for  them  could  be  derived. 
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i .  Introduction 

The  subject  known  as  viscoelasticity  has  developed  consider¬ 
ably  In  recent  years.  Th's  subject  deals  with  materials  whose 
stress-strain  behavior  is  time  dependent  for  all  levels  of  stress 
and  strain. 

The  mathematical  formulation  of  the  theory  of  viscoelasticity 
is  weU  established  (see  for  example  references  1-4).  The 
essence  of  any  such  formulation  is  the  appropriate  mathematical 
representation  of  the  stress-strain  behavior  of  the  viscoelastic 
material.  For  most  viscoelastic  materials  undergoing  small 
deformations,  the  dependence  of  stress  on  strain  can  be  well 
approximated  to  by  a  linear  relation.  In  this  case  there  are 
three  basic  approaches  which  are  commonly  used  for  describing 
the  stress-strain  behavior  of  a  viscoelastic  material. 

In  the  first  approach,  the  stress-strain  relation  is 
described  by  a  suitable  linear  differential  equation  with  constant 
coefficients  [2].*  In  order  to  describe  the  behavior  of  real 
viscoelastic  solids,  however,  such  an  equation  must  involve  quite 
a  few  terms  and  this  makes  any  mathematical  treatment  rather 
lengthy  and  in  fact  impractical.  In  the  second  approach,  the 
stress  strain  relation  is  described  by  a  suitable  hereditary 
integral  [1-4].  Such  a  representation  Is  based  on  the  Boltzmann 
superposition  principle  which  holds  for  linear  viscoelastic  solids. 
In  this  representation,  the  stress  is  evaluated  from  an  integral 
*  - 

Numbers  in  square  brackets  refer  to  the  bibliography  at  the 
end  of  the  paper. 
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which  includes  a  certain  material  property  known  as  the  stress 
relaxation  function  (or  relaxation  modulus),  and  the  strain  is 
evaluated  from  an  integral  which  includes  another  material 
property  known  as  the  creep  function  (or  creep  compliance).  The 
stress-relaxation  function  describes  the  stress-time  dependence 
of  the  material  when  subjected  to  a  unit  step  loading  of  strain. 
The  creep  function  in  contrast  describes  the  strain-time  depend¬ 
ence  of  the  material  when  subjected  to  a  unit  step  loading  of 
stress.  Both  properties  may  be  derived  for  any  linear  visco¬ 
elastic  material  from  fairly  simple  experiments  (see  for  example. 
Ferry,  pp .  105-150).  Such  experiments  may  usually  provide 
reliable  data  for  times  which  are  larger  than  1  sec.  For  times 
which  are  much  shorter  than  1  sec.,  the  stress-strain  behavior 
of  the  material  is  best  represented  by  the  third  approach. 

In  the  third  approach,  a  sinusoidal  stress  is  applied  to 
the  material,  and  the  resulting  sinusoidal  strain  is  related  to 
the  stress  through  a  complex  modulus  which  is  dependent  on  the 
frequency  of  the  applied  oscillations.  The  absolute  value  of 
this  complex  modulus  gives  the  ratio  between  the  maximum 
amplitudes  of  the  stress  and  the  strain,  whereas  the  ratio 
between  the  imaginary  and  real  parts  of  this  complex  modulus,  is 
a  measure  of  the  internal  friction  in  the  material.  There  are 
quite  a  few  experimental  methods  which  have  been  used  in  order  to 
determine  the  complex  modulus  of  viscoelastic  materials.  Two 
good  summaries  of  such  methods  were  given  by  Kolsky  [5-6]. 

In  this  report,  a  new  technique  is  described,  in  which 
the  complex  modulus  cf  viscoelastic  materials  (or  alternatively 
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the  phase  velocity  and  loss  angle  6  which  are  directly  related 
to  it  as  is  shown  in  appendix  I)  as  well  as  the  stress  relaxation 
function  E(t)  for  short  times,  may  be  evaluated.  This  technique 
is  based  on  certain  phenomena  which  take  place  during  the  arrival 
of  an  extensional  stress  pulse  at  the  interface  between  two  long 
bars  made  of  different  materials.  It  is  known  (see  for  example 
Kolsky  and  Lee  L  7  J )  that  an  extensional  stress  pulse  impinging 
on  the  boundary  between  two  elastic  rods,  will  split  into  two 
extensional  pulses,  namely  a  reflected  and  a  transmitted  pulse. 
The  amplitudes  of  the  reflected  and  transmitted  pulses  are 
related  to  that  of  the  incident  pulse  by  a  certain  characteristic 
impedance  ratio  which  is  evaluated  in  the  context  of  this  report 
(see  also  Kolsky  and  Lee  [7]).  However  the  shapes  of  all  three 
pulses  will  remain  similar. 

This  situation  no  longer  exists  when  one  or  both  of  the 
bars  are  viscoelastic.  In  this  case,  as  was  shown  by  Kolsky  and 
Lee  [7],  not  only  the  amplitudes  but  also  the  shapes  of  the 
reflected  and  transmitted  pulses  will  in  general  differ  from 
that  of  the  incident  pulse .  The  changes  in  pulse  shapes  are 
strictly  dependent  on  the  particular  properties  of  the  materials 
involved.  Hence  by  analysing  these  changes  by  a  suitable 
mathematical  procedure,  the  desired  viscoelastic  properties  may 
be  derived.  The  phenomenon  of  reflection  and  transmission  just 
described  and  its  implied  consequences,  form  the  basis  of  the 
experimental  technique  described  In  this  report. 


2 .  General  Description  of  Technique 

In  the  experimental  technique  described  in  this  report,  two 
bars,  one  of  which  was  elastic  with  known  properties  and  the 
other  viscoelastic  with  properties  to  be  determined,  were  placed 
end  to  end.  This  was  done  either  by  applying  a  thin  layer  of 
grease  between  the  two  bars  and  pressing  them  together  or  by 
gluing  one  bar  on  top  of  the  other. 

The  two  bars  were  held  in  a  vertical  position,  with  the 
elastic  bar  on  top  of  the  viscoelastic  bar.  A  small  steel  ball 
was  dropped  onto  the  top  surface  of  the  upper  bar.  The  impact 
produced  a  mechanical  pulse  which  was  propagated  along  this  bar 
towards  the  lower  bar.  At  the  interface  between  the  two  bars, 
reflected  and  transmitted  pulses  were  generated.  Recording  of 
all  three  pulses  (the  Incident  pulse,  the  reflected  pulse  and 
the  transmitted  pulse)  was  carried  out  by  means  of  electrical 
resistance  strain  gages  which  were  bonded  to  the  surface  of  the 
bars  at  different  locations.  The  strain  gage  outputs  were  fed 
on  to  a  cathode-ray  oscillograph  where  they  were  amplified  and 
the  traces  recorded  photographically.  The  photographic  images 
thus  obtained  were  first  magnified  by  the  use  of  a  slide  projector 
and  the  enlarged  pictures  were  traced  with  a  pencil. 

The  incident  and  reflected  pulse  shapes  were  later  analysed 
mathematically  with  the  aid  of  a  Fourier  analysis  or  alternatively 
by  a  Laplace  transform  method.  (It  Is  shown  in  this  report  that 
a  mathematical  analysis  of  the  transmitted  pulse  Is  not  required.) 
From  such  an  analysis,  as  is  shown  In  this  report,  the  unknown 


material  properties  of  the  viscoelastic  bar,  such  as  C  (the 
phase  velocity),  tan  5  (the  loss  factor)  and  E(t)  (the  stres 
relaxation  function)  could  theoretically  be  derived. 
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3 .  Theoretical  Background  and  Analysis 
3.1.  Elastic  Bars 

In  trying  to  find  material  properties  of  an  unknown 
material  by  means  of  the  experimental  technique  previously 
described,  let  us  start  our  investigation  with  the  relatively 
simple  case  of  two  elastic  bars  attached  together. 

Consider  an  extensional  wave  of  arbitrary  shape 
traveling  from  one  elastic  bar  to  another.  Upon  reaching  the 
boundary  between  the  two  bars,  this  wave  will,  in  general 
generate  a  reflected  pulse  and  a  transmitted  pulse.  Let  the  dis¬ 
placements  associated  with  the  incident,  the  reflected  and  the 
transmitted  pulses  be  denoted  by  u^,  u2  and  respectively. 

Assuming  the  incident  pulse  to  travel  in  the  direction 
of  increasing  x  where  x  is  the  longitudinal  coordinate  measured 
along  the  axis  of  the  bar,  we  may  then  write: 

u,  =  f(t  -  £-) 

1 

u  2  =  F(t  +  £-)  (1) 

u3  "  ®  ^  -  j=r^) 

where  and  are  the  wave  velocities  in  the  first  and  second 
bar  respectively.  The  relations  between  the  functions  f,  P,  and 
G  may  be  inferred  from  the  conditions  of  the  continuity  of  dis¬ 
placements  and  stresses  across  the  boundary  which  will  be  hither¬ 
to  designated  as  x=0  for  the  sake  of  convenience.  We  thus  have: 


U1  +  LU  -  U-, 


a  ^  +  a  2  =  c?^ 


at  x=0 


(Here  a  denotes  stress) . 

Applying  Hooke's  law  to  the  second  condition,  and  using  expression 
(1)  for  the  displacements  one  gets  at  x=0 

f(t)  +  F(t)  =  G(t) 

(2) 

El["  +  §YF’(t)]  =  "E2  ^,(t) 

where  the  dash  denotes  differentiation  with  respect  to  the 

argument.  Here  E-j^  and  E2  are  Young's  moduli  in  the  first  and 

2  El 

second  rod  respectively,  in  which  case  we  have  C,  ®  —  and 

2  E2  1 
C„  =  ,  p,  and  p0  being  the  corresponding  densities. 

P  2  i-  C- 

Differentiating  the  first  equation  in  (2)  w.r.t.  time  and  solving 
for  the  two  unknowns  F'(t)  and  G'(t),  one  gets: 


plCl”p2C2 

p  C  ^  t  p  2 


f  *  (t ) 


G'(t)  = 


2plCl 

P121+P2!22 


f '(t) 


Integrating  the  last  equations  w.r.t.  time,  and  assuming  zero 
initial  conditions,  one  gets: 


p  C ,  -  p  0C  p 

P(t)  f(t) 


2p  c 

G(t)  =  p~r"  +  o  p~ 
P1C l+p222 


f  (t) 
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Equations  (3)  give  us  the  existing  relations  between  the  three 
functions  f,  F,  and  G  which  describe  the  displacements  of  the 
pulses.  Observing  that  at  the  boundary  x=0,  =  f(t),  u2  =  F(t) 

and  u^  =  G(t),  we  may  deduce  the  following  relations  between  the 
displacement  amplitudes: 


u. 


u 


—  P2^2 
1  P1^1+P2C2 


2pi°i 


Pl^l+P2^2 


O) 


at  x=0  and  for  any  time  t  within  the  time  interval  at  which 
reflection  and  transmission  take  place  at  the  boundary.  In  the 
experimental  procedure,  however,  the  quantities  which  were  being 
measured  were  strains  rather  than  displacements. 

In  order  to  find  the  appropriate  relations  between  the 
strain  amplitudes,  let  us  first  derive  the  corresponding  strains 
from  expressions  (1)  and  (3).  We  then  will  have  at  x=0: 


3u^ 

3  x  x=o 


f  *  (t) 


£  2 


e-,  = 


d  u  2 

3x  j x=o 
3u-. 


F’(t)  = 


1_  P1C1~P2C2 
C1  P1C1+P2C2 


f  *  (t ) 


3  3  x  x  =  o 


G'(t)  =  - 


2p1C1 


C0  D,C,  +  p„C, 

C..  J.  1  Cl  C 


f  *  ( t  ) 


Thus  the  existing  relations  between  the  strain  amplitudes,  at 
x=0  and  for  all  times  within  the  time  1  erval  at  which  reflection 
and  transmission  take  place,  will  be 


lows  : 
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^2  P  9  —  p  2  ^  2  1 

e~  ~  ~  P1Ci+p2C2  “  P2C2+P1C1 


C,  2p1C1 

(  _±)  - ± — i 

^2  pl^l+p2^2 


(5) 


Expressions  (4)  and  (5)  show  us  that  the  various  relations  between 
the  pulses  are  determined  by  the  corresponding  characteristic 


impedances  pC  but  the  velocity  ratio  also  appears  in  the 

°2 

relations  for  strains. 


3 . 2 .  Viscoelastic  Bars 

Having  discussed  the  case  of  two  elastic  bars,  we  may 
now  proceed  to  the  more  complicated  case  of  two  viscoelastic  bars 
attached  together. 

As  a  starting  point  for  the  analysis  we  will  first 
consider  the  case  of  an  extensional  harmonic  wave  of  some  known 
frequency  traveling  from  one  viscoelastic  bar  to  another.  As  in 
the  case  of  the  two  elastic  bars,  this  situation  will  again  result 
in  the  generation  of  two  pulses  at  the  boundary,  namely  a 
reflected  pulse  and  a  transmitted  pulse.  Denoting  the  displace¬ 
ments  associated  with  the  incident,  the  reflected  and  the  trans¬ 
mitted  pulses  by  u^,  ur  and  u^_  respectively,  we  may  then  write 
(see  appendix  I )  : 

u1  =  A  exp[  -ot^x+ip  ( t  -  £r-)] 

ur  =  B  exp[+a1x  +  ip(t  +  ~)  ]  (6) 

ut  =  T  exp[  -otpX  +  ip  ( t  -  ~ )] 
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Here  and  are  the  attenuation  coefficients  in  the  first  and 
second  bar  respectively  and  and  C2  are  the  corresponding 
phase  relocities  at  the  given  frequency  p. 

A,  B  and  T  are  in  general  complex  so  that 

A  =  A1+iA2,  B  =  B1+1B2  ,  and  T  =  Th+iTg 

where  A1,A2,B1,B2,T^,T2  are  the  amplitudes  of  the  displacements 
at  the  interface  x=0. 

From  the  simple  one  dimensional  theory  of  wave  propaga¬ 
tion  in  viscoelastic  materials  (see  appendix  I),  we  have  the 
following  expressions  for  the  attentuation  a  and  the  phase 
velocity  C  in  a  viscoelastic  material: 


Here  p  is  the  density  and: 

tan  6  =  ~  E*  =  ^F^+E2 

where  E-L  +  1E9  is  the  complex  Young's  modulus  of  the  material  (see 
for  example  Ferry  [3]).  For  most  viscoelastic  materials,  commonly 
used  we  have 

tan'  6  <<  1 


In  this  case  we  may  write  the  following  approximations: 


a 


p  tan  6 
2C 


(7) 
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The  existing  relations  between  the  displacement  amplitudes 
A1,....,T2  may  again  be  found  by  applying  the  conditions  of  the 
continuity  cf  displacements  and  stresses  across  the  boundary. 
This  indeed  was  done  by  Kolsky  and  Lee  [73  who  found  the  follow¬ 
ing  relations  (see  Kolsky  and  Lee,  p.  25): 

!_s  2s(6  -6  ) 

B  =  ±— —  p  +  - ^ - L_  A 

1  1+s  1  (1+s)2  2 

2s ( 6 , -6  ~  )  ,  . 

B  =  - A,  +  k0 

2  ,.\i  1  1+S  2 


(l+s)' 

?  2s  ( 5?-6  ) 

T  =  — - —  4-  _ z _ i _ A 

1  1+5  (l+s)2  2 


(8) 


2s(<5  -6  ) 

T  =  - i - - —  A  +  — _  A 

2  ( !+s ) 2  1  1+5  2 


Here : 


s  = 


P2C2 

P1C1 


5,  = 


2  tan  6 


6.  = 


tan  6  ' 


where  <5  and  6'  are  the  phase  angles  between  a  harmonic  stress  and 
the  corresponding  harmonic  strain,  in  the  first  and  second  bar 
respectively . 

For  our  purposes,  however,  we  are  more  interested  in 
expressing  the  material  properties  of  the  second  bar  In  terms  of 
the  displacement  amplitudes,  since  the  latter  (as  will  be  shown 
later  in  the  analysis) are  directly  related  to  the  strains  which 
in  turn  will  be  measured  by  experiment. 


e 


i 
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In  order  to  find  an  expression  for  the  material  properties 
of  the  second  bar,  one  might  use  expressions  (8)  and  through 
some  algebraic  manipulations  arrive  at  some  complicated  formulas 
relating  6^  and  s  to  the  known  quantities  A1,A2,B1,B2,6^ . 

Another  way  of  finding  out  the  desired  relations  between  material 
properties  and  displacement  amplitudes  will  be  to  start  again  with 
the  same  analysis  p"^^^Fed  by  Kolsky  and  Lee  [7],  and  then  arrive 
directly  at  the  desired  relations  through  some  intermediate  step 
in  the  analysis.  The  second  approach  will  be  presented  here  as 
It  results  in  somewhat  simpler  relations  between  the  material 
properties  of  the  second  bar  and  the  displacement  amplitudes. 

We  may  now  start  our  analysis  by  writing  down  the  boundary 
conditions  for  the  displacements  and  stresses  across  the  boundary 
x=0 : 


UI  +  ur  =  ut 


fr  +  a  =  o. 
i  r  t 


at  x=0 


(9) 


The  stresses  o  a  and  o  may  be  derived  from  expression  (6) 

1  I  U 

using  the  following  relations: 

3u . 

°i  (E1  +  iE2) 

dur 

=  (El  +  iE2>  9^T 
°t  =  (E1+1E2>  3^ 


(10) 


Here  E^+iE0  Is  the  complex  Young's  modulus  In  the  first  rod  and 

»  t 

E.+iE^  Is  the  corresponding  modulus  In  the  second  rod.  Substitu¬ 
tion  of  the  displacements  and  stresses  from  expressions  (6)  and 


•%&*> 
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(10)  In  the  boundary  conditions  (9)  yields  the  following  relations 


(a)  A1  +  B1  =  Tx 

(b)  a2  +  b2  =  t2 


(ID 


(E1+iE2)[(-a1-  l£)(A1+iA2)  +  (a1+^-)  (B^+iBg)] 

’  <V1E2><-“2-  5f><TX+1T2> 


Now : 


E,  =  E^  tan  6 
E2  =  tan  6  • 


Also : 


p  tan  <5  ' 


2  tan  5 


2  tan  5  » 


EiCl 

C1C2 


P2C2 

P1C1 


If  we  now  make  these  substitutions  in  equation  (11c),  equate  real 

2 

and  imaginary  parts  and  neglect  terms  involving  tan  "6  w.r.t. 
unity,  we  obtain  the  following  two  relations: 


i4 


(A1~B1  )S  i  +  (A2_B2)  =  S(T162+T2) 
(A1-B1)  -  (A2~b2 ) 6 i  =  S(T1~T262) 


(12) 


Substitution  of  the  amplitudes  and  T2  from  expressions  (11a) 
and  (11b)  in  the  two  equations  (12)  yields: 


(A1-B1)61  +  (A2-B2)  =  S(A1+B1)62  +  S(A2+B2) 
(A1-B1)  -  (A2-B2)61  =  S(A1+B1)  -  S(A2+B2)62 


(13) 


The  last  two  equations  can  be  solved  algebraically  for  the  two 
unknowns  62  and  S.  When  this  is  done,  and  upon  replacing  6^  and 
6 2  by  i  tan  6  and  ^  tan  6'  one  gets: 


tan  S'  = 


S  = 


P2C2 

P1C1 


(A2+A2-B2-B2)tan  6  +  4  (A^-A^, ) 
(A2+A2-B2-B2)  +  (A1B2-A2B1)tan  6 

(A2+A2-B2-B2)  +  (A1B2-ApB1)tan  6 
(A1+B1)2  +  (A2+B2)2 


(14) 


Expressions  ( 14 )  give  us  the  desired  relations  between  the  material 
properties  of  the  second  bar,  namely  tan  6'  and  C0  and  the  dis¬ 
placement  amplitudes  of  the  incident  and  reflected  pulses. 

It  is  worth  noting  that  the  displacement  amplitudes  of 
the  transmitted  pu]se  do  not  enter  into  these  relations.  Hence 
In  carrying  out  experiments  which  are  directed  towards  the  evalua¬ 
tion  of  material  properties  of  the  second  bar,  ro  measurements  have 
to  be  taken  on  that  bar.  The  only  measurements  required  are 
those  obtained  from  the  first  bar-.  This  last  pact  Is.  of  extreme 
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Importance  when  one  is  trying  to  determine  mechanical  properties 
of  materials  of  very  large  damping  where  direct  propagation 
measurements  are  impractical,  or  materials  like  rubber  or  clay 
to  which  regular  measurement  techniques  (such  as  strain  gages) 
are  very  difficult  to  apply. 

Having  now  found  the  desired  relations  between  the 
material  properties  of  the  second  bar,  and  the  displacement 
amplitudes  and  material  properties  in  the  first  bar  for  a 
sinusoidal  wave  at  an  angular  frequency  p,  we  may  now  try  to 
evaluate  those  relations  for  the  case  of  a  general  pulse  of 
arbitrary  shape  and  duration.  In  order  to  do  this  we  must  first 
give  the  displacement  pulse  a  suitable  representation .  This 
might  be  done  by  writing  the  appropriate  pulse  in  terms  of  a 
Fourier  integral.  Thus  we  may  write  the  incident  pulse,  the 
reflected  pulse  and  the  transmitted  pulse  as: 

7  ,  ,  -a,x+ip(t-  — ) 

Uj  =  Re  (A,+iA„)e  1  dp 

JL  J  ±  d 

°  (1?) 

00  +a1x+ip(t+  ) 

u  =  Re  (B,+iB„)e  ^  dp 

r  i  d 

o 

00  -a  9  x  + 1  p  ( t.  - 

u t  =  Re  (T1  +  iT2)e  2  dp 

o 

The  displacements  described  by  expressions  (15)  are  no  more  than 
the  generalization  of  the  displacements  given  In  expressions  (6) 
for  the  case  of  an  infinite  range  of  continuous  frequency 
spectrum.  For  the  purpose  of  numerical  computations,  however. 

It  Is  necessary  to  represent  the  pulses  in  the  form  of  Fourier 


I 
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sums  rather  than  Fourier  integrals.  Thus  we  can  write  the 
three  pulses  as: 

•  -o1x+ikp0(t-  §-) 

u,  =  Re  E  (An+iAp)e 
1  k=o  L  ^ 

«  +alX+ikpo(t+  £-) 

u  =  Re  £  (B  +iB?)e  1  (16) 

r  k-o 

»  -<V+ikpo(t-  h~> 

u,  =  Re  I  (T,+iT~)e 
r  k=o  1  d 

In  the  last  representation  we  have  chosen  some  basic  frequency 

pQ  and  then  described  all  the  pulses  as  periodic  functions  with 
2  IT 

period  Tq  =  —  ,  and  Fourier  coefficients  A^,....,T2.  We  must 
therefore  observe  that  the  pulse  representation  given  by  ( 1 6 )  is 
not  the  same  as  that  given  by  (15).  The  latter  mathematically 
describes  a  pulse  whereas  the  first  one  describes  a  periodic 
function  with  period  T  . 

However,  within  the  time  Interval  from  sero  to  T  ,  both 
representations  are  equivalent,  namely  they  describe  exactly  one 
unique  pulse.  Hence,  In  using  expression  (16)  for  further 
analysis,  we  must  restrict  ourselves  to  times  which  are  not  less 
than  zero  and  not  greater  than  T  . 

The  basic  period  Tq  should  be  chosen  so  as  to  be  long 
enough  to  contain  the  whole  of  the  pulse,  but  also  small  enough 
to  insure  fast  convergence  of  the  Fourier  sum  to  its  final  value. 
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If  we  now  examine  the  displacements  described  by 

expressions  (16),  we  see  that  they  are  no  more  than  linear  sums 

of  expressions  similar  to  those  given  in  (6).  Since  we  are 

dealing  with  linear  viscoelasticity,  relations  (8)  and  (14)  must 

hold  for  each  set  of  displacement  amplitudes  (or  Fourier 

coefficients)  A^,....,^  corresponding  to  a  certain  discrete 

frequency  kp  .  Thus  in  order  to  find  the  material  properties  of 
o 

the  second  bar  in  the  case  of  two  general  incident  and  reflected 
pulses,  we  must  first  find  the  Fourier  coefficients  of  these  two 
pulses  by  a  Fourier  analysis  technique,  and  then  use  relations 
(14)  for  each  set  of  Fourier  coefficient  corresponding  to  a 
certain  discrete  frequency  kpQ.  In  this  manner  we  may  get  the 
material  properties  of  the  second  bar,  namely  tan  <5  *  and  as 
a  function  of  frequency  for  frequency  values  which  are  Integer 
multiples  of  the  basic  frequency  p  .  From  the  analysis  we  had 
so  far  we  could  find  the  existing  relations  between  the  material 
properties  of  the  second  bar  and  the  displacement  amplitudes  in 
the  first  bar.  However,  in  the  experimental  technique  described 
in  this  report,  the  quantities  being  measured  were  strains 
rather  than  displacements .  Hence  It  would  seem  appropriate  to 
carry  on  the  analysis  one  step  further  and  evaluate  the  relations 
between  the  material  properties  of  the  second  bar  and  the  strains 
in  the  first  bar. 

In  order  to  do  this  we  must  first  derive  appropriate 
expressions  for  the  strain  pulses  traveling  in  the  first  bar. 

This  might  be  done  by  referring  to  the  incident  and  reflected 


! 


displacement  pulses  described  in  expressions  (16)  and  differ¬ 
entiating  them  w.r.t.  x„ 

Denoting  by  and  the  resulting  incident  and 
reflected  strain  pulses,  we  may  then  write  at  x=0: 


Eih-o  =  Re  -  —  )(A1+lA2)e 


2  [  C-°iA1+ 

k=o  1 


g —  A2)cus  kp^t  +  (a^A^t  g -  A1)sin  kpQt] 


ev>L=~  -  Re  Z  (tor  +  —= - )(B,+iB0)e 


ikp^t 


(17) 


r  1  x^o 


1  C1  '  1  2 1 


00  kp  kp 

=  £  [  ( a^B1~  g—  B2)cos  kpQt  f  (-a1B2-  g—  B,)sin  ’Pct] 

k=o  1  u  '1 

Expressions  (17)  might  be  recognized  as  Fourier  series 

expansions  of  the  incident  and  reflected  strain  pulses  at  x=0. 

Thus  if  we  call  the  Fourier  cosine  coefficients  and  the 

Fourier-sine  coefficients  o;  the  incident  strain  pulse,  a  and 

K  c 

aks  respectively,  and  those  of  the  reflected  strain  pulse  bkc 
and  b.  respectively,  we  must  then  have: 

K  S 


-cm  A  .  +  ^ — -  A  _  -  a, 
lx  G,  2  kc 


aiA2  +  cy  ai  *  ak! 


(18) 


a1B1  -  jr--  B0  =  b. 
11  G  2  k  c 


-a  B0  -  ~ —  B -  b, 
12  C1  1  ks 
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Solving  the  first  pair  of  equations  for  the  displacement 
amplitudes  A1,A2,  and  the  second  pair  of  equations  for  the  dis¬ 
placement  amplitudes  B-^Bg,  and  also  using  the  relation 
(kp  )tan  6 

“i  *  — Jcl - one  gets: 

,  tan  6 

A  »  (!x.)  2  akc  _aks 

1  ^kp  ;  .  .  2 

0  !  +  (tan_6) 


_  .  tan  <5 

C,  a,  +  ^=5 —  a. 

/  1  v  kc _ 2 _ ks 

KP°  1  +  (^-g— )2 


c  tan-  6  b  _b 
fLl  a  2  kc  bks 

' kp  .  x  2 

°  i  +  ra-g-~) 


A  .  -KC  +  ^  V 

kP°  1+<^V 


(19) 


Inserting  the  last  expressions  for  the  displacement  amplitudes 
in  relations  (14)  and  doing  some  algebraic  manipulations,  we 
finally  get: 


(af  +af  -b,  -b2  )tan  6  -4(a,  b,  -a,  b,  ) 

.  ,,  kc  ks  kc  ks  kc  ks  ks  kc 

tan  o  '  =  — p - p - P - p - 

(akc+aks"bkc'bks)  +  (aw,bura^\Jtar>  6 


kc  ks  ks  kc 


P2C2 


(a.2  +af  -b,L  -b f  )  +  (a,  b,  -a,  b,  )tan  6  (£-0) 

kc  ks  kc  ks  K  im  v*  u *=  \,n 


^akc_bkc^  + 


kc  ks  ks  k c_ 

,  -b,  )2 

ks  ks 


Expressions  (20)  give  us  the  desired  relations  between  the 
material  properties  of  the  second  bar  and  the  strain  amplitudes 
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(or  Fourier  coefficients)  in  the  first  bar  for  all  frequencies 
kpo(k=l,2, . .  . ) . 


3.3.  Elastic  to  Viscoelastic  Bar 

In  the  theoretical  analysis  presented  so  far  we  have 
evaluated  some  basic  relations  which  exist  between  strain 
amplitudes  and  material  properties  in  one  bar,  and  material 
properties  of  a  second  bar  which  is  attached  to  the  first  one. 
When  both  bars  are  elastic,  these  relations  are  given  by  the 
first  expression  in  (5).  When  both  bars  are  viscoelastic,  these 
relations  are  given  by  expressions  (20).  When  the  first  bar 
(namely  the  one  in  which  the  incident  and  reflected  strain  pulses 
are  traveling)  is  elastic,  and  the  second  bar  is  viscoelastic, 
these  relations  can  be  obtained  from  expression  (20)  by  putting 
tan  6=0.  In  this  case  we  get: 


tan  o '  = 


"*  ^akcbks~aksbkc  ^ 
'  2  '  ?  2  2 
akc+aks~bkc"bks 


s  = 


P2C  2 
P1C1 


2.2  2  ,  2 

akc  aks  kc~  ks 


<akc‘bkc>2  +  (aks-bkl2 


(21) 


where  again  aRc  and  aRg  are  the  Fourier  cosine  coefficients  and 
Fourier  sine  coefficients  of  the  incident  strain  pulse,  and  b 

kc 

and  bks  are  tne  corresponding  coefficients  of  the  reflected 
strain  pulse.  When  the  first,  bar  Is  elastic  and  the  second  bai¬ 
ls  viscoelastic,  we  can  also  derive  an  additional  expression 


which  will  relate  the  material 


properties  and  strain  pulses  in 
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the  first  bar  to  another  material  property  of  the  second  bar, 
namely  the  stress  relaxation  function  E(t).  (For  the  definition 
of  E(t)  see  for  example  Flftgge,  p.  23.)  In  order  to  do  this  we 
must  solve  an  initial-boundary-value  problem  of  pulse  propagation 
in  a  system  of  two  bars  attached  together,  one  of  which  is 
elastic  and  the  other  viscoelastic. 

Consider  an  extensional  pulse  traveling  along  an  elastic 
bar  and  reaching  the  boundary  between  this  bar  and  a  visco¬ 
elastic  bar  attached  to  it.  As  was  already  mentioned  earlier, 
this  pulse  after  being  reflected  at  the  Interface  will  generate 
two  additional  pulses,  namely  a  reflected  pulse  and  a  trans¬ 
mitted  pulse.  Let-  us  again  denote  the  displacements  associated 
with  the  incident,  the  reflected  and  the  transmitted  pulses  by 
u  ,  u  ,  and  u  respectively.  The  equations  of  motion  for  all 
three  pulses  may  then  be  written  as: 


o 

9"u 


i 


3t“ 

~2 
3  u 


at 

„2 
a  u. 


at‘ 


i_  d_Z± 

P± 


1_  ^ 
Pj  9x 


1_  ^t 
3x 


(22) 


where  as  before,  o  denotes  stress,  p  denotes  density,  x 
denotes  distance  along  the  bar  axis,  and  t  denotes  time. 

For  the  sake  of  convenience  we  shall  consider  the  initial 
time  t=0  to  be  the  time  at  which  the  head  of  the  incident  pulse 
has  .just  reached  the  Interface,  and  the  two  other  pulses  have 
just  started  being  formed. 


The  positive  x  axis  will  be  chosen  to  point  in  the 
direction  of  the  viscoelastic  bar,  and  the  origin  x=0  will 
again  be  assumed  to  be  in  the  interface  between  the  two  bars. 
In  order  to  solve  the  equations  of  motion  (22)  we  must  first 
write  down  the  constitutive  relations  for  each  one  of  the 
three  pulses.  For  the  incident  and  reflected  pulses,  since 
they  are  propagating  in  an  elastic  medium,  these  relations  are 
simply  Hooke's  law,  namely: 


ai  E1  3x 


ar  =  E1  9x 


(23) 


Here  is  Young's  modulus  in  the  elastic  bar.  For  the  trans¬ 
mitted  pulse,  since  it  is  propagating  in  a  viscoelastic  medium. 


the  constitutive  relation  may  be  formulated  in  terms  of  a 
hereditary  integral  (see  for  example  Fldgge,  p.  24),  namely: 


u  9e  u  9  u 

at(t)  =  {  E2(t-T)  5^  dr  -  J  E2(t-r)  ^  dr  (2H) 


Here  et  is  the  transmitted  strain  and  E2(t-x)  Is  the  stress 

relaxation  function  in  the  viscoelastic  bar  at  the  time  t-x. 

If  all  the  pulses  with  which  we  are  dealing  are  assumed  to  be 

9  c. 

smooth  functions  of  time,  then  ^ -  Is  zero  for  all  tmes 

9  x 

between  and  0+  and  hence  relation  (24)  may  be  written  as: 


OtU)  =  E2(t-T)  Jjpf  dT 


(25) 


,v* '  • ;  ‘  m 
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Inserting  now  the  constitutive  relations  (23)  and  (25)  in  the 
equations  of  motion  (22)  one  gets: 

E. 


2 

8  u . 


3t2  P1  3  x2 


2  2 
1  a  ui  ^  ui 

=  C1  T“2“ 


3x 


3  2u 


3 1* 


E1  „2  a2ur 

P-,  „  2  C1  .  2 


i  3  x 


3x 


(26) 


„  2 
3  u, 


3t‘ 


3 _ 

3x 


32u, 

E2(t~T)  3x3t 


dr  = 


t 

u 

nH 


3^u. 


E2(t-T) 


3x23t 


dx 


Here  is  the  pulse  velocity  in  the  elastic  bar  which  is  known 
to  be  equal  to  /E^/p^  Equations  (26)  will  now  be  solved  by 
means  of  a  Laplace  transform  technique.  Taking  the  Laplace 
transform  of  the  second  equation  In  (26)  and  denoting  by  p  the 
transform  variable,  and  by  the  transformed  reflected  dis- 

x 

placement,  one  gets: 


P 


or : 


(27) 


Taking  now  the  Laplace  transform  of  the  third  equation  in  (26) 
and  denoting  by  the  transformed  transmitted  displacement  and 
by  E2  the  transformed  stress  relaxation  function,  one  gets : 


P 


3 


2- 

u 


h 


3  x 


2 


inannim 


or: 
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„2— 

3  u.  p0 

— —  pu,  =  0  (28) 

3x2  E2  c 


In  the  derivation  of  equations  (27)  and  (28),  zero  initial 
condition  were  assumed,  which  is  in  accordance  with  the  assump¬ 
tion  already  stated  that  all  pulses  are  considered  to  be  smooth 
functions  of  time. 

Equations  (27)  and  (28)  can  now  be  solved  for  the  dis¬ 
placements  u  and  u.  thus  yielding: 
r  T/ 


Assuming  the  two  bars  to  be  very  long,  we  must  have: 


The  condition  =  0  Is  necessary  in  order  to  prevent  an 
infinite  growth  of  the  reflected  pulse  on  its  traveJ  In  the 
negative  x  direction.  The  condition  B-,  =  0  is  necessary  in 

-L 

order  to  prevent  an  infinite  growth  of  the  transmitted  pulse  on 
its  travel  in  the  positive  x  direction.  Thus,  expressions  (29) 


may  be  reduced  to: 
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(30) 


The  two  constants  and  B ^  might  be  found  from  the  boundary 
conditions  at  x=Q.  Those  conditions  require  the  equality  of 
displacements  and  stresses  across  the  boundary,  and  therefore 
may  be  written  as: 


(a)  u1  +  Uj,  -  u1 


(b)  +  ar  «  at 


at.  x=0 


(31) 


Here  ar  and  cr^  denote  transformed  stresses.  Expressing  now 

the  constitutive  relations  (23)  and  (25)  in  terms  of  their 
Laplace  transforms  we  get: 


3u 

o  =  E,  r — 
r  1  3x 


_  _  3u, 

°t  =  E  ? 


(32) 


Inserting  the  last  expressions  for  the  transformed  stresses  in 
the  second  boundary  condition  (31b;  one  gets: 


9Uj 


l'3x 


3  u 

_ i 

3  x 


-)  -  E0P 


3u . 


3  x 


at 


x  =  0 


(31c) 


Rewriting  now  the  boundary  conditions  (31a)  and  (31c)  In  a 
slightly  different  form  we  get: 
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u  =  u 
r  i 


at  x=Q 


(33) 


j  P  3  U| 

P  3x“ 


3u_  3u. 


3x  3x 


Now,  since  the  incident  pulse  is  propagating  in  the  elastic  bar, 
we  know  that  it  has  to  be  expressed  in  the  form: 


u±(x,t)  =  u1(t  -  ~) 


Taking  the  Laplace  transform  of  both  sides  we  get: 


x 

-  "  ^P- 
u4  =  e  u, (p) 

-L  J- 


(34) 


where  u^Cp)  is  the  transform  of  the  displacement  shape  u^(t). 

Differentiating  expression  (34)  w.r.t.  x,  we  also  get: 


3Uj 

air 


_  E_  e 

ci 


Vp) 


(35) 


Before  using  the  boundary  conditions  (33)  we  also  need  to  find 

3u  3u, 

expressions  for  ^ — -  and  ~ —  . 

o  X  o  X 

These  may  be  found  by  differentiating  w.r.t  x  the  dis¬ 
placements  in  expression  (30),  thus  yielding: 


3  u 


3x 


r  ^  P_  A  °1 


C1  "1 


(36a) 


:  Pr—ri » 


! 

I 

i 


27 


3ut  _ 

n - 


B2e 


(36b) 


Now  we  can  substitute  expressions  (30),  (3*0,  (35)  and  (36)  in 
the  boundary  conditions  (33),  thus  getting  (at  x-0)  : 


B 


2  ~  A1  =  Uj 


_?2p  g 

"E1  P/e2 


(37) 


*>  b2  -  §7  Ai 


-  ui(p) 


Solving  the  last  two  equations  for  the  constants  and  ,  and 

2 

substituting  p^C1  for  E,,  one  gets: 


A,  = 


Bo  * 


)1C1"  /P2E2  P 


1/2 


u,  (p) 


P1C1+  y/p2E2  pl/2 

u,  (p) 


2PiCi 


(38) 


P1E1+  ^2^2  p 


Inserting  the  constants  and  from  expressions  (38)  in 
expression  (30)  for  the  transformed  displacements,  one  gets 


(a) 


-  C7X  Pl°l  Jp2E2  P 

u  =  e  1 - 


l/'< 


>lCl+  >/b  2E 


u  1/2 


u1(p) 


J¥px 


(39) 


(b) 


ut  =  6 


2plCl 


p  1C1+  Vp7e9  p‘ 


n1  (p) 
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Expressions  (39)  give  us  the  existing  relations  between  the 
reflected  and  transmitted  pulses  and  the  incident  pulse,  in 
terms  of  transformed  quantities . 

Observing  now  that  the  reflected  pulse  is  traveling  in  the 
elastic  bar,  it  must  have  the  following  representation: 

ur(x,t)  =  ur(t  +  JJ-) 

Taking  the  Laplace  transform  of  both  sides  we  get: 


x 

-  - 

ur  =  e  ur(p) 


(40) 


where  ur(p)  is  the  transform  of  the  displacement  shape  u^(t). 

Inserting  the  last  expression  for  ur  in  the  left-hand 
side  of  expression  (39a),  we  get: 


1  — 
e  ur(p)  =  e 


Cxx  p^-  /pgE 


i?  .1/2 


3ici+  pl/2 


u±(p) 


or: 


ur(p)  P1C1 


-/ 


P2E2  p 


1/2 


u,  (p) 


-1/2 


(41) 


pici+  y p2e2  p 


Expression  (4l)  gives  us  the  relation  between  the 
material  properties  in  the  elastic  and  viscoelastic  bars  anl 
the  transformed  displacement  shapes  in  the  elastic  bar.  How¬ 
ever,  from  an  experimental  point  of  view  we  are  more  interested 
in  finding  a  relation  similar  to  ( )  which  will  involve  strain 
shapes  rather  than  displacement  shapes.  In  order  to  do  this  let 


>v\* «  ,  ... 
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us  write  again  the  expressions  for  the  incident  and  reflected 
displacement  pulses  in  the  elastic  bar: 

U±(x,t )  -  u1(t  -  ~) 
ur(x,t)  =  ur(t  +  £-) 


Differentiating  the  last  expressions  once  w.r.t.  x  and  once 
w.r.t.  t,  we  get: 

3Ui  _  _  1_  '  Xk 

3x  ei  “  ~  ui  ' v 


=  vj  =  u{(t  -  §-) 


3u, 

9F 

au. 


3x  er 


k  u-(t  +  k] 


(42) 


3u* 

at 


=  vr  =  ur(t  +  g-) 


Here  e^, > er > vr  are  the  strains  and  particle  velocities  in  the 
incident  and  reflected  pulse... 

From  expressions  (42)  the  following  relations  are 
obtained : 


3Uj 

W 


-  V± 


3u 

r 

at  ' 


Vr 


(43) 


Relations  (43)  will  hold  for  any  location  x  in  the  bar.  If  we 

now  fix  our  attention  to  some  arbitrary  fixed  point  X  in  the 

3u.  3u 

bar,  then  the  quantities  ^ which  appear  in 

expressions  (43)  may  be  all  considered  as  functions  of  time  only 
at  that  fixed  point  X. 
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Hence  we  can  integrate  expressions  (43)  w.r.t.  to  time 
and  get: 


u,  (t)  =  -C-,  [  eA 

±  1  j  - L 


t »  )dt 


t 

ur(t)  =  C-,  j  er(t '  )dt 


(44) 


Taking  now  the  Laplace  transform  of  both  sides  of  each  one  of 
the  expressions  in  (44)  we  get: 


uj  (p) 


ur(p) 


er(p) 

P 


(45) 


Here  e^(p)  and  er(p)  denote  transformed  strain  shapes  correspond¬ 
ing  to  e^(t)  and  e  (t).  Inserting  now  the  last  expressions  for 
u^p)  and  ur(p)  in  expression  ( 4 1 ) ,  we  get  * 

e  (p)  p,  C,  -  / p„E~  pP//l  </p„E„  p1/2-PlC. 

zr - =  -  — ^f=± - 7  =  — — - —  (46) 

ei(p)  P1C1+  /p^2  P1/2  7 P2E2  p1/2+p1C] 

Expression  (46)  gives  us  the  desired  relations  between  the 
material  properties  In  the  elastic  and  viscoelastic  bars  and  the 
transformed  strain  shapes  In  the  elastic  bar.  It  is  worth  noting 
that  when  the  viscoelastic  bar  is  replaced  by  an  elastic  bar, 
the  stress  relaxation  function  E9(t')  turns  out  to  be  simply  the 
constant  Young's  modulus  E0.  In  this  case: 
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e2  =  r 


Inserting  this  value  for  E2  in  expression  (46)  we  get: 

er(p)  2^2  ”pl^l  p2^2~Pl^l 
e^Cp)  +p1C1  p2C2+plCl 

Taking  the  inverse  Laplace  transform  of  the  last  expression  we 


get : 


£p(t)  (>2^2  pl^l 

(t)  “  P2c2  +  pici 


The  last  expression  for 


e„(t ) 


is  the  same  as  the  first  relation 


in  (5)  which  was  obtained  directly  for  the  case  of  two  elastic 
bars  attached  together.  Returning  to  the  problem  of  a  visco¬ 
elastic  bar  attached  to  an  elastic  bar,  we  would  like  to  have 
an  explicit  expression  for  E2  in  terms  of  all  the  other  related 
quantities  which  appear  in  our  problem.  This  might  be  done  with 
the  aid  of  expression  (46),  yielding: 


( p  C, )2  e.+e 

b2  =  <=^> 

t-i  *-r 


(4?) 


If  we  now  take  the  Inverse  Laplace  transform  of  E0  which  was 
found  by  using  expression  (47),  the  stress  relaxation  function 
E2(t)  can  be  evaluated. 

The  necessary  steps  In  the  evaluation  of  E„(t)  are  there¬ 


fore  as  follows  : 


First,  the  incident  and  reflected  strain  pulses  e^t)  and  er(t) 
in  the  elastic  bar  are  recorded  graphically.  Next,  the  trans¬ 
formed  strains  Fj(p)  and  if  (p)  are  computed  by  using  the  regular 
Laplace  transform  formulas,  namely: 

CO  'p 

e^(p)  =  j  e_ptei I't )dt  =  e-ptei(t)dt 
o  o 

T  ^8> 
er,(p)  =  e“pter(t)dt  =  j  e“pter(t)dt 
o  o 

Here  T  is  the  duration  of  the  two  puls.s,  which  accounts  for 
the  change  in  the  upper  limit  of  integration  from  00  to  T.  Since 
the  strains  e ^ ( o )  and  er(t)  are  not  known  in  an  analytic  form, 
expressions  (48)  have  to  be  evaluated  numerically.  This  could 
be  done  for  example  by  evaluating  the  integrals  in  (48)  accord¬ 
ing  to  Simpson’s  rule  using  different  values  for  p.  In  this 
way  the  values  of  e^(p)  and  er(p)  could  be  obtained  for  any 
chosen  value  of  p.  The  next  step  will  be  to  substitute  the 
transformed  strains  ^(p)  and  e”r(p)  in  expression  (47),  thus 
getting  the  value  of  E'2  for  any  chosen  value  of  p.  The  last 
step  and  certainly  the  most  difficult  one  is  to  find  E^Ct)  as 
an  inverse  Laplace  transform  of  E0 .  We  must  observe  in  this 
respect,  that  finding  the  inverse  Laplace  transform  of  E2 
would  not  be  an  easy  task  even  if  E.,  would  have  been  given  in 
the  form  of  an  analytic  expression  in  p .  In  our  case  the 
situation  is  even  more  complicated  owing  to  the  fact  that  is 
given  in  the  form  of  point  values  at  different  data  poines  p, 
rather  than  in  the  form  of  an  analytic  function  in  the  variable  j 


Fortunately  however ,  for  most  known  viscoelastic  materials, 

Eg(t)  does  not  change  very  rapidly  w.r.t.  time  (see,  for 
example.  Ferry,  Fig.  2-2  in  p.  29).  In  such  a  case  we  can  find 
the  inverse  Laplace  transform  by  a  very  nice  and  simple  approxi¬ 
mate  method  developed  by  Schapery  [8],  According  to  this  method 
E2(t)  can  be  found  from  Eg  by  the  following  approximate 
formula: 

Eg ( t )  =  [pEg]  Q  R  (49) 

p=  -~- 

That  is.  Eg ( t )  3s  found  by  multiplying  the  corresponding  trans¬ 
formed  function  by  the  transform  variable  and  evaluating  the 

0  5 

result  at  points  p  corresponding  to  points  . 

Since  the  function  Eg  is  given  only  for  discrete  data 
points  p,  also  the  result  E^t)  will  be  calculated  only  for 
discrete  time  values  t.  Thus  by  finding  Eg  from  the  experi¬ 
mental  values  of  e^(t)  and  er(t)  using  expressions  (48)  and  (47), 
and  by  later  applying  to  it  the  approximate  inversion  formula 
(49),  the  stress  relaxation  function  Eg(t)  may  be  evaluated  for 
discrete  time  values  t. 
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4 .  Experimental 

4 . 1  Experimental  Set-up 

The  experimental  set-up  which  was  used  for  the  experi¬ 
ments  described  in  this  report  is  shown  in  Figs,  la,  lb  and  1c . 

The  oscilloscope  shown  in  this  set-up  was  a  Tectronix 
dual-beam  scope  type  502A,  and  it  had  a  polaroid  camera  attached 
to  its  screen.  The  triggering  of  the  scope  was  done  by  the 
dropping  ball,  which  upon  hitting  the  upper  bar,  closed  an 
electrical  triggering  circuit.  The  triggering  circuit  used  is 
shown  in  Fig.  2. 

When  the  upper  bar  was  of  a  conducting  material,  it 
formed  part  of  the  triggering  circuit.  When  the  upper  bar  was 
an  insulator,  a  thin  layer  of  conducting  paste  was  applied  to 
the  top  face  of  the  bar,  thus  enabling  the  dropping  ball  to  close 
the  triggering  circuit  upon  hitting  the  edge  of  the  bar. 

The  bars  which  were  used  in  the  experiments  were  in 
general  2  feet  long  and  1/2"  in  diameter.  They  were  usually 
resting  on  the  table,  but  sometimes  suspended  with  the  aid  of 
3  strings  coming  down  from  the  ceiling.  The  dropping  ball  which 
was  made  of  steel,  was  generally  1/2"  in  diameter. 

The  strain  gages  used  in  the  experlmenis  were  paper  base 
wire  gages  of  the  type  SR-4,  C-10  made  by  the  BLH  Company. 

These  gages  had  a  resistance  of  1000Q  and  gage  factor  cf  3.0  and 
this  made  them  capable  of  detecting  very  low  strains  of  the  order 
of  magnitude  of  10  .  Each  of  the  strain  gages  which  was  used 

for  measuring  the  transient  strain  pulses,  was  connected  in  a 
dynamic  bridge  circuit  as  shown  In  Fig.  3. 
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The  input  voltage  E  to  this  circuit,  was  produced  by  a 
set  of  6V  batteries  joined  in  series.  When  a  strain  pulse 
passed  through  the  measuring  gage,  the  output  voltage  E  was 

fr 

changed  by  the  amount  dE  which  was  proportional  to  the  strain 

8 

e.  The  voltage  change  dE  was  fed  to  the  scope  and  recorded 

o 

photographically  by  the  camera.  It  is  worth  noting  that 
conversion  of  the  output  voltage-change  dE  to  strain  units, 

o 

although  possible,  was  not  in  general  necessary.  The  reason 
for  this  is  that  in  the  formulas  which  relate  the  material 
properties  to  the  strain  amplitudes,  only  strain  ratios  and 
not  absolute  strain  values  are  Involved.  Since  the  strains  are 
proportional  to  the  corresponding  output  voltage-changes  dE  , 

o 

their  ratio  is  the  same  as  the  ratio  of  the  corresponding 
voltage-changes.  The  measuring  stiain  gages  were  bonded  to  the 
bars  using  either  Eastman  910  adhesive  or  Duco  cement.  The 
strain  gage  In  the  upper  bar  was  bonded  at  a  distance  of  about 
9^-"  from  the  Interface,  and  was  thus  capable  of  giving  separate 
recordings  of  the  incident  and  reflected  pulses  without  any 
Interference  between  them.  The  upper  bar  was  always  chosen  to 
be  elastic  or  very  nearly  so,  and  therefore  there  was  no 
appreciable  difference  between  strain  pulses  being  recorded  at 
the  strain  gage  location,  and  the  same  pulses  while  traveling 
near  the  Interface.  This  fact  is  important  since  in  the 
formulas  relating  and  tan  6 '  to  the  Fourier  coefficients  of 
the  incident  and  reflected  strain  pulses,  all  coefficients 
involved  are  those  belonging  to  strain  pulses  acting  at  the 
boundary  x  =  0 .  Also,  In  the  analysis  which  relates  E.,(t )  to  the 
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incident  and  reflected  strain  shapes,  it  was  assumed  from  the 
beginning  that  the  first  bar  was  elastic.  The  other  bar  which 
could  be  either  elastic  or  viscoelastic,  had  at  least  one  strain 
gage  mounted  on  it  for  recording  the  transmitted  strain  pulse. 
This  recorded  pulse,  though  not  required  for  direct  computations 
of  material  properties  in  the  lower  bar  was  used  for  additional 
information  and  as  a  check  on  the  results.  When  the  second  bar 
was  elastic,  it  was  sufficient  to  glue  to  it  one  strain  gage  at 
any  distance  from  the  interface.  Since  no  dispersion  took  place 
in  this  case,  the  recorded  pulse  had  the  same  shape  as  it  had 
while  traveling  near  the  interface  or  at  any  point  in  the  bar. 
When  the  lower  bar  was  viscoelastic,  two  strain  gages  were 
employed;  one  was  mounted  as  close  as  possible  to  the  interface, 
and  the  other  was  mounted  at  some  distance  from  the  interface . 

The  first  strain  gage  recorded  the  transmitted  strain  pulse  near 
the  Interface,  immediately  after  it  had  been  generated,  and  the 
second  strain  gage  recorded  the  same  pulse  after  it  had  travelled 
some  distance  in  this  bar  and  had  undergone  considerable 
dispersion . 

The  dummy  strain  gages  which  were  used  in  the  dynamic 
bridge  circuits,  were  bonded  to  a  piece  of  metal  and  put  in  a 
closed  metal  box  which  provided  the  necessary  shielding  against 
external  electrical  noise.  The  attachment  of  the  two  bars  was 
achieved  either  by  applying  a  thin  layer  of  grease  between  them 
and  pressing  them  together  using  a  short  close-fit  metal  collar 
to  insure  their  stability,  or  by  gluing  them  together  with 
Eastman  910  cement.  Before  such  attachment  was  made,  the  two 
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edges  to  be  connected  were  first  machined  on  a  lathe  and  then 
thoroughly  lapped  using  a  specially  designed  lapping  tool. 

This  lapping  tool  is  shown  diagramatically  in  Fig.  4. 

The  lapping  operation  was  aone  by  putting  some  compound 
paste  between  the  metal  block  and  the  steel  cylinder  and  rubbing 
the  one  against  the  other.  In  some  cases  it  was  found  that  a 
slight  roughening  of  the  bar  edges  was  necessary  in  order  to  get 
better  adhesion  between  the  bars.  In  such  cases.  Instead  of 
compound  paste,  very  fine  emery  paper  was  attached  to  the  metal 
block,  which  in  turn  was  rubbed  against  the  steel  cylinder  as 
before.  All  the  electrical  wires  in  the  experimental  set-up 
(excluding  the  very  short  ones  only)  were  shielded  wires;  this 
was  done  in  order  to  prevent  external  electrical  noise  from 
interfering  with  the  outputs  of  the  strain  gages. 

4 . 2  Experimental  Procedure 

The  bars  which  were  used  in  the  experiments  were  made 
of  the  following  materials:  Steel,  aluminum,  polystyrene , 
polyethylene,  Lucite  (p.m.m),  green  solid  propellant  (AVCO 
product),  pinewood  (direction  of  grain  parallel  to  axis  of  bar), 
sandstone,  compressed  sand,  and  clay.  (The  densities  of  these 
materials  are  given  In  appendix  II.) 

Before  the  experimental  procedure  (which  will  later  be 
described)  was  used  with  these  bars,  some  preliminary  experi¬ 
ments  were  conducted  the  object  of  which  was  to  check  various 
aspects  of  the  overall  experimental  technique;  these  experiments 
will  be  described  first. 
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The  first  preliminary  experiment  was  to  hit  a  single 
steel  bar  with  a  dropping  ball  and  record  the  output  of  the 
strain-gage  which  was  attached  to  it.  The  result  of  this  experi¬ 
ment  is  shown  in  Fig.  5-  The  object  of  this  experiment  was  to 
check  the  functioning  of  the  overall  experimental  set-up,  and  to 
see  what  kind  of  a  strain  pulse  was  generated  in  the  bar  upon 
its  being  hit  by  the  dropping  ball. 

The  second  preliminary  experiment  was  done  by  joining 
together  two  similar  steel  bars,  with  the  aid  of  grease  and  a 
short  metal  collar,  and  hitting  the  upper  bar  with  a  dropping 
ball.  The  outputs  of  the  strain  gages  attached  to  both  bars 
were  recorded,  and  are  shown  in  Fig.  6.  The  object  of  this 
experiment  was  to  check  whether  or  not  good  contact  between  two 
attached  bars  could  be  achieved  by  this  means. 

The  next  preliminary  experiments  were  done  by  joining 
together  an  aluminum  and  a  steel  bar  with  a  grease  joint,  and 
performing  the  same  experiment,  first  with  the  aluminum  bar  on 
top,  and  then  with  the  steel  bar  on  top.  The  results  of  these 
experiments  are  shown  in  Figs.  7  and  8.  The  object  of  these 
experiments  was  to  check  the  agreement  between  theory  and  experi¬ 
ment  in  the  case  of  materials  with  known  properties,  and  thus 
provide  an  additional  verification  of  the  adequacy  of  the  experi¬ 
mental  technique,  to  evaluate  material  properties  of  new 
materials . 

home  additional  preliminary  experiments  were  then  carried 
out  with  polystyrene  bars,  since  it.  was  decided  to  use  a  poly¬ 
styrene  bar  as  the  upper  bar  in  most,  of  the  experiments  described 
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in  this  report.  This  choice  was  made  since  polystyrene  is  very 
close  to  being  elastic,  and  its  characteristic  impedance  does 
not  differ  greatly  from  the  characteristic  Impedances  of  most  of 
the  other  materials  whose  properties  were  being  sought.  This 
last  consideration  is  important  since  the  ratio  between  the 
characteristic  impedances  of  the  upper  and  lower  bars  must  be 
such  that  while  the  amplitude  of  the  reflected  pulse  should  not 
be  too  small  compared  with  that  of  the  incident  pulse,  it  should 
also  not  be  too  close  to  it  in  magnitude. 

Thus  if  the  reflected  strain  pulse  be  too  small,  it  will 
be  difficult  to  distinguish  it  from  the  small  disturbances 
which  are  present  in  the  tail  of  the  incident  pulse.  If  on  the 
other  hand,  the  reflected  pulse  has  almost  the  same  amplitude 
(in  absolute  value)  as  the  incident  pulse,  large  errors  may 
appear  in  the  computations  of  the  material  properties  of  the 
lower  bar. 

The  reason  for  this  is  that  these  computations  are 
largely  dependent  on  the  arithmetic  difference  between  the 
absolute  values  of  the  incident  and  reflected  strain  amplitudes. 
This  could  be  shown  for  example  by  referring  to  expressions  (5) 
which  apply  to  the  relatively  simple  case  of  two  elastic  bars 
attached  together.  from  the  first  expression  in  (5)  we  get: 
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Now,  if  both  strains  have  the  same  sign,  trie  difference?  between 
their1  absolute  values  appears,  In  the  denominator  of  the  last 


expression.  If  on  the  other  hand,  one  strain  has  the  opposite 
sign  of  the  other,  the  difference  between  their  absolute  values 
will  appear  in  the  numerator  of  the  last  expression.  At  any 
rate,  we  see  that  the  material  property  of  the  second  bar, 
namely  C2,  depends  largely  on  the  difference  between  the 
absolute  values  of  and  er •  If  this  difference  is  too  small, 
any  slight  change  in  one  of  the  strains  will  considerably  effect 
the  difference  between  the  strains,  and  hence  cause  enormous 
errors  in  the  computation  of  C2.  In  some  of  the  experiments 
described  in  this  report,  it  was  indeed  necessary  to  choose  for 
upper  bars,  steel  or  aluminum  bars  rather  than  polystyrene  bars  . 
This  was  done  in  order  to  get  a  better  match  between  the 
characteristic  impedance  of  the  upper  bar  and  that  of  the  lower 
bar,  according  to  the  considerations  discussed  above. 

Another  reason  for  choosing  most  of  the  upper  bars  to  be 
made  of  polystyrene  rather  than  another  glass-like  plastic  is, 
that  a  considerable  amount  of  data  is  available  on  this  material. 
The  1  irst  experiment  on  polystyrene  was  to  hit,  a  single  poly¬ 
styrene  bar  with  the  dropping  ball,  and  the  result  obtained  is 
shown  in  Fig.  9.  The  objects  of  this  experiment  were:  1)  to 
see  what  kind  of  a  pulse  is  produced  in  these  circumstances, 

2)  to  determine  the  velocity  of  propagation  of  such  a  pulse  in 
the  bar,  3)  to  check  how  close  the  behavior  of  polystyrene 
approximated  to  perfect  elasticity  since,  as  mentioned  above, 
it  was  necessary  for  the  upper  bar  to  be  elastic  or  very  close 
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The  second  experiment  on  polystyrene  was  done  with  two 
polystyrene  bars  bonded  with  Eastman  910  cement .  The  top  bar 
was  hit  by  a  dropping  ball;  the  result  of  this  experiment  Is 
shown  in  Fig.  10.  The  object  of  this  experiment  was  to  check 
whether  or  not  it  was  possible  to  use  cement  as  a  means  for 
attachment  of  the  two  bars  without  introducing  appreciable 
errors  in  the  experimental  results. 

Having  described  so  far  the  preliminary  experiments 
which  were  conducted  on  some  of  the  bars,  we  will  now  describe 
the  main  experimental  procedure  which  was  applied  to  different 
bars  in  order  to  evaluate  their  material  properties. 

The  first  step  in  the  procedure  was  to  put  an  elastic 
or  nearly  elastic  bar  on  top  of  the  bar  whose  properties  were 
being  sought.  When  the  lower  bar  was  made  of  po’ yethylene ,  green 
solid  propellant,  wood,  compressed  sand,  or  clay,  the  upper  bar 
was  chosen  to  be  a  polystyrene  bar.  When  the  lower  bar  was  made 
of  Lucite,  the  upper  bar  was  chosen  to  be  an  aluminum  bar. 

When  the  lower  bar  was  made  of  sandstone,  the  upper  bar  was 
chosen  to  be  a  steel  bar. 

The  particular  choice  of  the  elastic  bar  in  each  case, 
was  such  that  its  characteristic  impedance  matched  that  of  the 
lower  bar,  according  to  the  considerations  mentioned  earlier 
in  connection  with  the  choice  of  polystyrene  .  The  next  step  in 
the  experimental  procedure  was  to  hit  the  upper  bar  with  the 
dropping  ball  a  number  of  times,  and  take  the  recordings  from 
the  different  strain  gages  mounted  on  the  bars,  for  different 
settings  of  amplification  and  beam  sweep  rate  in  the  scope. 
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The  results  obtained  in  this  way  for  different  kinds  of 
materials  are  all  presented  in  Pigs.  11  to  32.  The  last  step 
in  the  experimental  procedure  was  to  take  one  of  the  pictures 
showing  the  incident  and  reflected  strain  pulses,  project  it  on 
a  sheet  of  paper  attached  to  the  wall  and  trace  the  magnified 
picture  with  a  pencil.  The  magnified  figures  thus  obtained 
could  be  analyzed  according  to  the  theory  presented  in  the 
previous  chapter.  The  results  of  such  an  analysis  should  have 
yielded  the  desired  properties  of  the  lower  bar.  However,  as 
will  be  shown  in  the  next  chapter,  the  information  which  was 
obtained  in  this  way  was  very  limited  owing  to  the  large 
sensitivity  of  the  calculated  results  to  small  changes  in  the 
pulse  shape  and  the  pulse  starting  time.  Such  changes  may 
always  be  present  as  a  result  of  small  inaccuracies  which  are 
inherent  in  any  experimental  procedure. 

4 . 3  Experimental  Results 

We  may  start  with  the  result,  of  the  first  preliminary 
experiment  which  is  the  picture  shown  in  Fig.  5.  This  picture, 
which  was  obtained  by  hitting  a  single  steel  bar  with  the 
dropping  ball,  shows  the  same  strain  pulse  travelling  back  and 
forth  in  the  bar,  and  changing  from  compression  to  tension  or 
vice  versa,  each  time  reflection  takes  place  at  one  of  the  free 
edges  of  the  bar.  Similar  pictures  were  obtained  whether  the 
bar  was  resting  on  the  table,  suspended  vertically  in  the  air, 
or  suspended  horizontally  in  the  air  and  hit  by  a  bail  coining 
from  the  side  . 
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Th~  duration  cf  the  strain  pulse  shown  in  Fig.  5  Is 
about  30  psec.  The  second  preliminary  experiment  which  was  done 
on  two  steel  bars  attached  together,  gave  the  picture  shown  in 
Fig.  6.  In  thi  ;  picture  the  strain  pulse  in  the  upper  beam  is 
the  incident  strain,  and  the  first  pulse  in  the  lower  beam  is 
the  transmitted  strain.  The  other  pulses  in  the  lower  beam  are 
the  same  transmitted  strain  being  intermittently  reflected  from 
the  bottom  edge  and  the  top  edge  of  the  lower  bar. 

Nov/,  since  the  two  bars  in  this  experiment  were  made  of 
the  same  material,  if  the  contact  was  perfect  it  would  be 
expected  that  the  incident  stress  pulse  which  was  travelling  in 
the  upper  bar,  would  be  entirely  transmitted  to  the  lower  bar 
and  no  reflection  should  occur  at  the  interface;  that  this 
indeed  proved  to  be  the  case  can  be  seen  from  Fig.  6. 

From  this  figure  we  also  note  that  the  grease  layer 
between  the  bars  acted  as  a  rectifier  for  stress  pulses  .  It 
allowed  the  passage  of  a  compression  pulse  from  the  upper  bar  to 
the  lower  bar,  but  prevented  the  passage  of  a  tension  pulse  from 
the  lower  bar  back  to  the  upper  bar.  This  of  course  is  due  to 
the  ability  of  grease  to  transmit  compression  but  yield  In 
tension,  causing  the  two  bars  to  separate  upon  the  arrival  of  a 
tension  pulse  at  the  interface.  Such  a  tension  pulse  in  the 
lower  bar  will  therefore  be  reflected  at  the  interface  as  a 
compression  pulse  in  that  bar,  since  the  Interface  at  that  time 
acts  as  a  free  boundary.  The  results  of  the  next  preliminary 
experiments  which  were  perforat'd  on  an  aluminum  bar  attached  to 
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a  steel  bar,  are  shown  in  Figs.  7  and  8.  Figure  7  was  obtained 
with  the  aluminum  bar  on  top,  and  Fig.  8  was  obtained  with  the 
steel  bar  on  top. 

In  both  figures  the  incident  and  reflected  pulses  are 
shown  in  the  upper  beam  and  the  transmitted  pulse  is  shown  in 
the  lower  beam.  If  the  maximum  amplitudes  of  those  pulses  are 
measured  directly  on  the  pictures,  it  is  found  that  the  relations 
between  them  (including  the  appropriate  sign,  which  can  be 
chosen  as  (+)  for  compression  and  (-)  for  tension,  or  vice  versa) 
are  in  very  good  agreement  with  the  theoretical  relations  given 
by  expressions  (5).  This  of  course  provides  an  additional 
verification  of  the  adequacy  of  the  experimental  technique  in 
evaluating  material  properties  from  strain  measurements,  for 
materials  with  unknown  properties. 

1 

In  Figures  7  and  8  we  again  see  that  after  the  reflected 
and  transmitted  strain'  pulses  have  been  generated,  the  Interface 
starts  to  act  as  a  free  boundary  in  each  one  of  the  bars,  owing 
to  the  Inability  of  the  grease  layer  to  sustain  tension. 

The-  next  preliminary  experiments  which  were  performed  on 
polystyrene  bars,  gave  the  results  shown  In  Figs.  9  and  10. 

Figure  9  which  was  obtained  by  hitting  a  single  polystyrene  bar 
with  the  dropping  ball  is  very  similar  to  Fig.  5  which  was 
obtained  by  performing  the  same  experiment  on  a  single  steel  bar. 
However  there  are  some  differences  between  the  results  of  these 
two  experiments.  In  the  polystyrene  bar  experiment,  the  duration 
of  the  strain  pulse  is  seen  to  be  about  180  usee  whereas  1i  is. 
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only  30  ysec  in  the  steel  bar  experiment.  Also  the  strain 
amplitudes  obtained  in  the  polystyrene  bar  experiment,  are  very 
much  higher  than  those  obtained  in  the  steel  bar  experiment  . 

These  differences  are  to  be  expected  from  the  differences  in 
Young's  modulus  between  polystyrene  and  steel  which  result  in 
corresponding  differences  in  both  duration  and  strain  as  given 
bv  hertz's  theory  of  indentation. 

Another  result  which  may  be  obtained  from  Pig.  9  is  the 
pulse  velocity  in  the  polystyrene  bar.  This  velocity  can  be 
calculated  from  the  time  it  took  the  pu1 ~e  to  travel  a  fixed 
distance  in  the  bar.  By  taking  several  measurements  on  Fig.  9> 
which  correspond  to  different  distances  of  travel,  the  velocity 
range  in  polystyrene  was  found  to  be  6000-6120  ft/sec.  For 
practical  purposes,  a  mean  velocity  of  6060  ft/sec.  or  ±850  m/sec. 
could  be  taken  as  the  pulse  velocity  in  polystyrene,  that  is: 

C  ,  .  =  6060  ft/sec.  +  1%  -  i860  m/sec.  +  1% . 

polystyrene  —  — 

Another  result  which  may  be  derived  from  Fig.  9  is  that  the 
attenuation  of  the  pulse  in  the  polystyrene  bar  is  indeed  very 
small.  Hence,  if  a  polystyrene  bar  is  used  as  the  upper  bar  in 
experiments  which  are  done  on  materials  like  polyethylene  or 
p.m.m.,  it  can  be  considered  to  be  elastic  for  all  practical 
purposes.  This  was  later  verified  by  some  additional  computa¬ 
tions  which  were  done  in  connection  with  the  derivation  of 
material  properties  ol  polyethylene .  These  computation:.;  will  be 
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The  result  of  the  second  experiment  with  polystyrene, 
namely  the  one  which  was  done  on  two  polystyrene  bars  bonded 
together,  is  the  picture  shown  in  Pig.  10.  Prom  this  picture 
we  see  that  good  contact  between  the  two  bars  had  been  achieved, 
and  the  whole  incident  pulse  was  transmitted  to  the  lower  bar, 
just  as  in  the  case  of  the  two  attached  steel  bars.  However, 
there  is  some  difference  between  the  two  cases .  In  the  steel 
bars  experiment,  the  transmitted  pulse  could  not  travel  back  to 
the  upper  bar  because  the  grease  layer  in  the  interface  could 
not  transmit  tension  pulses  (this  was  explained  before  in 
detail).  However  in  the  polystyrene  bars  experiment,  the  trans¬ 
mitted  pulse  could  travel  back  to  the  upper  bar  as  can  be  seen 
from  Fig.  10.  The  reason  is  that  the  cement  in  the  interface 
could  sustain  tension  as  well  as  compression,  and  hence  could 
transmit  the  tension  pulse  in  the  lower  bar  back  to  the  upper 
bar. 

Another  conclusion  which  may  be  derived  from  the  last 
experiment  is,  that  the  cement  probably  does  not  appear  to 
interfere  very  much  with  the  experimental  results,  and  can 
therefore  be  safely  used  as  a  means  for  the  attachment  of  one 
bar  to  another.  Having  presented  so  far  the  results  of  the 
preliminary  experiments,  we  will  now  present  the  results  of  the 
main  experiments,  which  for  the  sake  of  convenience  will  be 
described  for  each  material  separately. 

P  o.lyet:  hy  lone 

The  results  of  a  series  of  experiments  in  which  poly¬ 


styrene  bars  were  attached  to  polyethylene  bars ,  are  shown  in 
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Figs.  11,  12a,  12b,  and  13.  Figure  11  shows  the  incident  and 
reflected  strain  pulses  in  the  polystyrene  bar  (shown  in  the 
upper  beam)  and  the  transmitted  strain  pulse  in  the  polyethylene 
bar  (shown  in  the  lower  beam).  Figures  12a  and  12b,  which  were 
obtained  from  two  similar  experiments,  give  magnified  pictures 
of  the  incident  and  reflected  strain  pulses.  Figure  13  shows 
the  transmitted  strain  pulse  in  the  polyethylene  bar;  this  pulse 
was  recorded  twice  by  the  same  strain  gage,  once  while  traveling 
away  from  the  interface,  and  later  while  traveling  back  towards 
the  interface  after  being  reflected  from  the  bottom  edge  of  the 
bar . 

Figures  12a  and  12b  must  now  be  analyzed  using  the  theory 
presented  in  chapter  3  in  an  attempt  to  evaluate  the  material 
properties  of  polyethylene.  This  was  done  by  the  use  of  computer 
programs,  which  included;  1)  a  sub-program  for  the  evaluation 
of  the  Fourier  coefficients  of  the  incident  and  reflected  strain 
pulses,  2)  a  sub-program  for  tne  evaluation  of  the  Laplace 
transforms  of  the  incident  and  reflected  strain  pulses. 

3)  Formulas  for  the  evaluation  of  the  desired  material  properties, 
that  is,  expressions  (20) were  used  for  the  evaluation  of  tan  6' 
and  COJ  and  expressions  (97)  and  (99)  for  the  evaluation  of  E,,(t). 

c. 

In  using  these  programs,  a  basic  interval  of  2'JQ  gsec 
(corresponding  to  a  basic  frequency  of  3-7  kc/sec)  was  chosen 
for  each  one  of  the  pulses  in  Figs.  12a  and  12b.  Each  basic 
interval  was  divided  into  7b  divisions,  and  the  corresponding 
pulse  ordinates  were  read  into  the  programs  which  w<cre  later  run 
->n  the  computer.  It  turned  out  however,  that  most  of  the  results 
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obtained  in  this  way  were  rather  unreliable.  In  fact,  when  the 
ordinates  of  Fig.  12a  were  read  into  the  computer  program, 
tan  6*  was  obtained  negative  for  axl  the  Fourier  harmonics. 

When  the  ordinates  of  Fig.  12b  were  put  into  the  computer  tan  <5  ' 

was  obtained  positive  for  some  Fourier  harmonics  (including  the 

1st  harmonic),  and  negative  for  others,  and  it  always  differed 
very  much  from  one  harmonic  to  another.  All  these  results  are 
of  course  unacceptable,  since  it  is  known  that  tan  6*  for 
polyehtylene  is  about  0.1  and  is  more  or  less  constant  over  a 
wide  range  of  frequencies.  It  was  therefore  suspected  that  small 
inaccuracies  in  the  pulse  shape  and  also  in  the  location  of  the 
coordinate  system  used  for  the  calculation  of  the  Fourier 
coefficients,  could  cause  enormous  errors  in  the  results  of 
tan  6'.  This  indeed  proved  to  be  the  case  as  will  be  shown  in 
the  next  chapter.  It  is  worth  mentioning  that  the  results  of 

tan  6',  ridiculous  as  they  were,  were  not  affected  much  by  small 

variations  in  tan  6  (6  is  the  phase  angle  in  the  polystyrene  bar). 

In  fact  they  were  almost  the  same  for  tan  5  between  zero  and 
0.01  (the  correct  value  of  tan  6  is  less  than  0.01).  This 
verifies  the  assumption  that  polystyrene  could  be  considered  to 
be  elastic  for  all  practical  purposes. 

The  results  which  were  obtained  for  Ch,  were  quite  reason¬ 
able  for  the  first  harmonic  (the  basic  frequency),  but  again 
unreliable  for  the  higher  harmonics,  the  reasons,  being  the  same 
as  those  mentioned  above  in  connect i on  with  the  results  of  tan  2'. 

The  values  of  C.,  which  were  obtained  for  the  basic 
frequency  from  Figs .  12a  and  12b  were  920  m/ser  and  900  m/sec 

i 
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respectively.  (The  values  of  C2  which  were  obtained  for  the 
higher  harmonics,  were  sometimes  much  above  and  sometimes  much 
below  these  values.  This  is  unacceptable  since  according  to  the 
linear  theory  of  viscoelasticity,  a  linear  increase  of  C2  with 
the  logarithm  of  the  frequency  would  be  expected  when  tan  6  '  is 
constant,  and  this  is  more  or  less  so  for  polyetylene  . ) 

It  is  interesting  to  compare  these  values  with  the  value 
of  0  which  was  htained  by  Hillier  [9]  from  direct  experiments 
with  polyethylene.  Hillier’s  results  were  obtained  for 
sinusoidal  waves  of  different  frequencies  and  at  different  room 
temperatures  (see  Ref.9>  table  1,  p.  705).  The  relevant  part  of 
Hillier 's  results  is  again  presented  in  Appendix  III  for  the 
convenience  of  the  reader. 

Interpolating  Hillier's  results  for  a  frequency  of 
3-7  kc/sec  (which  was  the  basic  frequency  in  all  the  previous 
calculations)  and  a  temperature  of  13°C  (which  was  the  tempera¬ 
ture  in  tne  room  when  the  experiments  on  polyethylene  were 
performed),  the  resulting  value  for  C-,  turns  out  to  be  880  m/sec. 
This  result  is  in  good  agreement  with  the  results  of 
previously  obtained. 

Another  way  of  obtaining  an  approximate  value  for  C0  for 
the  first  harmonic  from  our  experiments  is  to  measure  the 
distance  between  the  pulse  peaks  in  Pig.  13  and  find  the  'pulse 
velocity'  5n  the  same  way  as  was.  done  for  the  polystyrene  bar. 

The  location  of  the  strain  gage  which  gave  the  result  shown  In 
Pig.  13,  was  bjj"  from  She  bottom  edge  of  the  pol  y- ■  1  hy  lone  bar. 
hi  nee  t,ho  distance  between  the  pulse  peaks  is  '■> .  0  cm  on  the 
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picture,  and  the  sweep  rate  is  0.1  millisec/cm,  the  ’pulse 

velocity*  turns  out  to  be  870  m/sec.  Now,  in  the  case  of  an 

elastic  bar’,  the  value  of  the  'pulse  velocity  '  is  close  to  the 

value  of  the  group  velocity  of  the  dominant  frequency  of  its 

dC 

Fourier  components.  This  group  velocity  is  equal  to  C-A  , 

where  C  is  the  phase  velocity  and  A  is  the  wave  length.  When 
dC 

A  -rj  is  much  smaller  than  C,  the  value  of  the  group  velocity  is 
almost  equal  to  the  value  of  the  phase  velocity.  Such  a  condi¬ 
tion  exists  for  example  when  the  pulse  length  is  much  greater 
than  the  diameter  of  the  bar,  in  which  case  the  dispersion  due 
to  geometrical  effects  can  be  neglected.  In  the  case  of  a 
viscoelastic  bar,  the  above  arguments  are  somewhat  harder  to 
justify  since  the  concept  of  group  velocity  is  not  well  defined. 
Nevertheless  it  seems  reasonable  tc  assume  chat  in  this  case 
also,  the  'pulse  velocity'  is  still  close  to  the  phase  velocity 
of  the  dominant  Fourier  harmonic,  when  A  ~  is  much  smaller 
than  C.  For  such  a  condition  to  exist,  (tan  6) /it  must  be  small 
compared  with  unity  (the  proof  of  this  statement  is  given  in 
Appendix  IV).  The  geometrical  dispersion  must  also  be  small  as 
In  the  elastic  case  and  this  applies  when  the  pulse  length  Is 
much  greater  than  the  diameter  of  the  bar.  The  conditions  just 
mentioned  were  indeed  valid  for  the  pulse  which  was  propagated 
in  the  polyethylene  bar. 

lienee,  the  value  which  was  found  for  the  'pulse  relocity’ 
(8/0  m/sec)  is  probably  very  close  to  the  value  of  the  phase 
velocity  C for  the  first  harmonic.  Thus  i  f  the  value  of  C,  for 

^  C- 

the  first  harmonic  is  considered  to  be  870  m/sec  to  the  first. 
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approximation,  we  see  that  this  value  is  in  very  good  agreement 
with  the  values  of  C2  previously  obtained  from  the  main 
theoretical  analysis  and  from  Hillier's  results.  Still  another 
way  of  finding  an  approximate  value  for  C0  for  the  first  harmonic 
from  our  experiments,  is  to  measure  the  maximum  strain  amplitudes 
in  Fig.  12a  or  Fig.  12b  and  use  expression  (50)  for  the  evalua¬ 
tion  of  Cg. 

Now,  expression  (50)  was  derived  for  two  elastic  bars 
attached  together.  However,  as  will  be  shown  in  the  next  chapter, 
the  shape  and  amplitude  of  the  reflected  pulse  will  be  almost 
the  same  whether  the  second  bar  is  elastic  or  viscoelastic  with 
small  tan  6'.  Hence  the  value  of  u2  derived  by  using  expression 
(50)  and  assuming  that  the  two  bars  are  elastic,  is  probably 
very  close  to  the  real  value  of  C0  for  the  first  harmonic  in  the 
viscoelastic  bar. 

Evaluating  C 2  for  the  first  harmonic  in  this  way  from 
Fig.  12b  (which  is  believed  to  be  the  more  reliable  figure 
between  the  two  figures  12a  and  12b),  it  is  found  that  its  value 
is  900  m/sec.  This  value  of  is  in  very  good  agreement  with 
all  the  other  values  of  previously  obtained  by  different 
methods . 

The  results  which  were  obtained  for  the  stress  relaxation 
function  E^(  ),  by  analysing  Figs.  12a  and  12b  according  to  the 
theory  presented  in  the  previous  chapter,  were  rather  unreliable, 
;;gain  due  to  the  large  effect  of  small  experimenta]  inaccuracies 
on  the  computed  results.  In  fact,  when  Fig.  Ira  was  used  in 
order  to  calculate  the  values  of  E^Ct),  for  different  values  of 
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time  t,  a  monotonically  increasing  function  w.r.t.  time  was 
obtained,  which  of  course  is  physically  unrealistic . 

When  small  shifts  were  made  in  the  location  of  the 
coordinate  system  of  the  reflected  pulse  (the  coordinate  system 
was  used  in  order  to  calculate  the  Laplace  transform  of  the 
pulse),  different  values  for  E2(t)  were  obtained,  but  the 
tendency  to  increase  with  time  was  still  unchanged.  When 
Pig.  12b  (which  has  slightly  different  pulse  shapes  than  those 
of  Fig.  12a)  was  used  in  the  calculations,  E-2(t)  was  obtained 
as  a  slightly  decreasing  function  w.r.t.  time.  However,  E0(t) 
could  become  an  increasing  function  w.r.t  time  even  in  this 
case,  if  a  small  shift  in  the  location  of  the  coordinate  system 
of  the  reflected  pulse  was  made.  Thus  we  see  that  the  results 
for  E2(t)  were  neither  reproducible  nor  reliable.  However  the 
average  value  of  E2(t)  was  very  close  to  some  effective  dynamic 
Young's  modulus  E^^,  whi  ch  can  be  set  equal  to  p2C2  where  C2 
can  be  obtained  by  any  one  of  the  methods  mentioned  above.  If 
the  approximate  value  of  C2  is  taken  as  900  m/sec,  the  apprcxi- 

O  P 

mate  value  of  E^^,  for  polyethylene  turns  out  to  be  7.9X10  kg/cm". 

A  mathematical  proof  showing  the  large  effects  of  small  changes 

In  coordinate  location  and  pulse  shapes  on  the  function  E0(t) 

£_ 

is  presented  in  the  next  chapter. 

Prom  the  results  obtained  for  polyethylene,  we  conclude 
that  the  only  quantity  which  can  be  determined  from  the  experi¬ 
ments  with  reasonable  accuracy,  is  the  phase  velocity  C2  at  the 
basic  frequency  p  in  the  viscoelastic  bar.  However  it  was 
shown,  that  in  order  to  find  this  phase  velocity,  there  was  no 
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need  to  use  the  rigorous  analysis  and  go  to  the  computer. 

Instead,  direct  measurements  on  Figs.  1.2  (either  12a  or  12b) 
and  13,  and  straightforward  calculations,  result  in  a  value  of 
C g  which  is  not  less  accurate  than  the  value  of  C2,  which  Is 
obtained  by  the  use  of  rigorous  analysis.  Thus,  in  order  to 
find  the  phase  velocity  for  materials  other  than  polyethylene, 
only  simple  methods  and  straightforward  calculations  will  be 
employed . 

Lucite  (p  .in ,m .  ) 

The  results  of  a  series  of  experiments  which  were  performed 

on  aluminum  bars  attached  to  Lucite  bars,  are  shown  in  Figs.  14, 

15  and  16.  Figure  14  shows  the  incident  and  reflected  strain 

pulses  in  the  aluminum  bar  (lower  beam)  and  the  transmitted 

strain  pulse  in  the  Lucite  bar  (upper  beam).  Figure  15  shows  a 

magnified  picture  of  the  incident  and  reflected  strain  pulses  in 

the  aluminum  bar.  Figure  16  shows  the  transmitted  pulse 

recorded  by  two  strain  gages  bonded  on  the  Lucite  bar  at  a 

distance  of  9^-"  from  each  other.  (The  upper  beam  corresponds  tc 

the  strain  gage  which  was  closer  to  the  interface.)  In  order  to 

find  an  approximate  value  for  the  phase  velocity  C0  in  the 

Lucite  bar,  the  strain  amplitudes  in  Fig.  15  were  measured  and 

inserted  in  expression  (50).  The  value  of  C2  thus  obtained  was 

2150  m/sec.  Another  approximate  value  for  was  calculated 

from  the  time  it  took  the  transmitted  puise,  shown  In  Fig.  16, 

to  travel  the  distance  between  the  two  recording  strain  gages  . 

The  value  of  C0  obtained  in  this  way  was  2070  m/spo. 

<- 
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The  value  of  the  correct  phase  velocity  at  the  basic 
frequency  p  ,  Is  probably  not  very  far  from  either  one  of  the 
two  values  previously  obtained,  and  hence  could  be  taken  approxi¬ 
mately  as  2100  m/sec.  The  value  of  the  effective  dynamic 
Young’s  modulus,  based  on  the  last  value  for  Cg,  turns  out  to 
be  5.3*10^  kg/cm^. 

Green  Solid  Propellant  (g.s.p.) 

The  results  of  a  series  of  experiments  done  on  a  poly¬ 
styrene  bar  attached  to  a  g.s.p.  bar  are  shown  in  Figs.  17,  18 
and  19.  Figure  17  shows  the  incident  and  reflected  pulses  in 
the  polystyrene  bar  (lower  beam)  and  the  transmitted  pulse  in 
the  g.s.p.  bar  (upper  beam). 

Figure  18  shows  a  magnified  picture  of  the  incident  and 
reflected  pulses  in  the  polystyrene  bar.  Figure  19  shows  the 
transmitted  pulse  twice  recorded  by  the  same  strain  gage  which 
was  located  at  a  distance  of  11”  from  the  bottom  edge  of  the 
g.s.p.  bar. 

The  value  of  C 2  which  was  found  from  Fig.  18  by  measuring 
the  strain  amplitudes  and  using  expression  (50),  turned  out  to 
be  1080  m/sec . 

The  value  of  which  was  found  from  Fig.  19,  by  dividing 

the  distance  of  travel  of  the  transmitted  pulse  by  the  time  of 

traverse,  turned  out  to  be  1050  m/sec.  If  the  correct  phase 

velocity  for  the  first  harmonic  is  assumed  to  have  the  value 

1060  m/sec.,  the  corresponding  effective  Young's  modulus  turns 

i|  2 

out  to  be  1.12*10  kg/cm  . 
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Pine  Wood 

The  results  of  a  series  of  experiments  which  were  carried 
out  with  a  polystyrene  bar  attached  to  a  pine  wood  bar  are 
shown  in  Pigs.  20,  21,  and  22.  The  result  of  one  experiment 
which  was  done  by  hitting  a  single  pine  wood  bar  with  the 
dropping  ball,  is  shown  in  Pig..  23. 

Figure  20  shows  the  incident  and  reflected  pulses  in  the 
polystyrene  bar  (upper  beam),  and  the  transmitted  pulse  in  the 
wooden  bar  (lower  beam).  Figure  21  shows  a  magnified  picture 
of  the  incident  and  reflected  strain  pulses.  Figure  22  shows 
the  transmitted  pulse  twice  recorded  by  a  strain  gage  which  was 
located  at  a  distance  of  15^"  from  the  bottom  edge  of  the 
wooden  bar.  Figure  23  shows  one  pulse  traveling  back  and  forth 
in  a  single  pine  wood  bar.  The  values  of  C2  which  were  found 
from  Figs.  21  and  22  in  the  same  way  as  for  the  g.s.p.  bar, 
were  5550  m/sec  and  5350  m/sec  respectively.  The  value  of 
which  was  found  f^nm  Fig.  23  in  the  way  as  was  done  in  tha 

single  polystyrene  bar  was  5500  m/sec. 

From  Figs.  22  and  23  we  see  that  the  decay  in  the  pulses 
traveling  in  the  wooden  bars  was  very  small,  and  hence,  wood  can 
be  considered  as  an  elastic  material  for  all  practical  purposes. 
If  we  assume  that  the  correct  value  of  the  wave  velocity  in  wood 
is  5500  m/scc..,  we  get  for  Young's  modulus  (by  writing  E  =  pC  ) 
the  value  2.16*10  kg/cm  which  is  in  good  agreement  with  the 
value  of  Young's  modulus  given  in  the  literature  for  that  kind 
of  wood  along  the  grain. 


Sandstone 
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The  results  of  a  series  of  experiments  with  a  steel  bar 
attached  to  a  sandstone  bar,  are  shown  in  Pigs.  2 4 ,  25  and  26. 

Figure  24  shows  the  incident  and  reflected  pulses  in 
the  steel  bar  (upper  beam)  and  the  transmitted  pulse  in  the 
sandstone  bar  (lower  beam).  Figure  25  shows  a  magnified  picture 
of  the  incident  and  reflected  pulses  in  the  steel  bar.  Figure 
26  shows  the  transmitted  pulse  twice  recorded  by  a  strain  gage 
which  was  located  at  a  distance  of  4^-g-  from  the  bottom  edge  of 
the  sandstone  bar. 

The  values  of  the  phase  velocity  C2  found  from  Figs.  25 

and  26  by  the  same  methods  as  before,  are  both  4100  m/sec. 

From  Fig.  26  it  can  be  seen  that  the  decay  in  the  sandstone  is 

very  small,  and  hence  sandstone  can  be  considered  to  be  elastic 

for  all  practical  purposes.  The  value  of  Young's  modulus  for 

5  2 

sandstone  turns  out  to  be  4.63*10  kg/cm  . 

Clay 

The  results  of  a  series  of  experiments  performed  on  a 
polystyrene  bar  attached  to  a  clay  bar  are  shown  in  Figs .  27, 

28  and  29.  Figure  27  shows  the  incident  and  reflected  strain 
pulses  in  the  polystyrene  bar  (upper  beam)  and  the  transmitted 
pulse  in  the  clay  bar  (lower  beam).  Figure  28  shows  a 
magnified  picture  of  the  incident  and  reflected  strain  pulses . 
Figure  29  shows  the  transmitted  pulse  recorded  by  two  strain 
gages  attached  to  the  clay  bar.  These  strain  gages  were  put 
8  cm  apart,  and  the  distance  between  the  lower  gage  and  the 
bottom  end  of  the  clay  bar  was  10  cm.  The  upper  beam  in  Fig.  29 
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corresponds  to  the  upper  gage,  and  the  lower  beam  corresponds 
to  the  lower  gage.  The  strain  gages  which  were  put  on  the  clay 
bar,  were  pressed  against  the  clay  arid  held  In  position  with  a 
small  piece  of  scotch  tape.  The  outputs  of  these  gages  are 
probably  far  from  giving  the  correct  values  of  the  strains  in 
the  clay,  since  the  strain  gages  are  much  stiffen  than  the  clay 
and  therefore  are  not  stretched  to  the  same  amount  as  the  clay. 
Nevertheless,  some  useful  information  can  still  be  obtained  from 
the  outputs  of  these  gages,  as  will  be  shown  rater. 

From  Fig.  28,  using  expression  (50),  the  phase  velocity 
in  clay  is  found  to  be  460  m/sec.  Another  value  of  C2  is  found 
by  referring  to  Fig.  29  and  dividing  the  distance  between  the 
two  recording  gages  by  the  time  it  took  the  transmitted  pulse  to 
travel  that  distance.  The  value  of  found  in  this  way  is 
4(0  m/sec. 

The  two  values  which  were  obtained  for  the  wave  velocity 
in  clay  are  in  good  agreement  with  the  value  obtained  by 
H.  Calvit  [10]  which  was  500  m/sec.  The  effective  dynamic 
Young's  modulus  based  on  a  phase  velocity  of  470  m/sec.  turns 
out  to  be  4.05*10^  kg/cm2. 

From  Fig.  29  it  can  be  seen  that  the  decay  of  a  strain 
pulse  in  clay  is  much  larger  than  it  is  in  the  plastic  materials 
previously  mentioned.  This  conclusion  is  mainly  based  on  the 
lower  beam  record  which  shows  the  transmitted  pulse  twice 
recorded  by  the  same  strain  gage.  Therefore,  even  if  the  pulses 
in  the  lower  beam  do  not  represent  the  true  strains,  the  ratio 
between  them  is  probably  close  to  the  ratio  between  the  true 
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strains  since  they  are  both  measured  by  the  same  strain  gage. 
Compressed  Sand 

The  experiments  with  sand  bars  were  done  on  wet  sand 
bars  and  on  dry  sand  bars.  In  both  cases,  the  bars  were  pre¬ 
pared  by  pressing  wet  sand  into  hard  paper  tubes  and  later 
unwrapping  these  tubes.  Some  of  the  sand  bars  which  were 
produced  in  this  way,  were  put  in  a  heated  oven  until  they 
became  dry. 

The  densities  of  the  wet  b^.rs  and  the  dry  bars  were 
found  to  be  about  the  same.  However  the  wet  sand  bars  were 
much  softer  than  the  dry  bars,  and  this  had  a  considerable 
effect  on  the  experimental  results  as  shown  below. 

The  results  of  a  series  of  experiments  done  on  poly¬ 
styrene  bars  attached  to  wet  sand  bars  are  shown  in  Pig.  30. 

This  figure  shows  the  incident  and  reflected  strain  pulses  in 
the  polystyrene  bar.  By  measuring  the  strain  amplitudes  on 
this  figure  and  using  expression  (50),  the  wave  velocity  in  wet 
sand  is  found  to  be  240  m/sec.  The  corresponding  Young's 
modulus  is  found  to  be  1.23*10^  kg/cm2. 

The  results  of  a  series  of  experiments  done  on  polystyrene 
bars  attached  to  dry  sand  bars,  are  shown  in  Figs.  31  and  32. 

Figure  31  shows  the  incident  and  reflected  strain  pulses  in  the 
polystyrene  bar.  Figure  32  shows  the  transmitted  pulse  in  the 
sand  bar  being  recorded  by  two  strain  gages  which  were  care¬ 
fully  bonded  to  the  bar  at  a  distance  of  3.5  cm  apart. 

From  Fig.  31,  and  by  using  expression  (50),  the  wave 
velocity  in  dry  sand  was  found  to  be  4 7 0  m/sec.  From  Fig.  32, 
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by  measuring  the  distance  between  the  two  pulse  peaks,  the  wave 
velocity  in  dry  sand  was  found  to  he  500  m/sec.  (The  strain 
gage  outputs  in  Pig.  32  do  not  represent  the  true  strains  since 
the  strain  gages  are  much  stiffer  than  the  sand.  However  this 
fact  does  not  affect  the  result  which  was  obtained  for  the  wave 
velocity.)  The  value  of  Young's  modulus  in  dry  sand  correspond¬ 
ing  to  an  average  wave  velocity  of  *180  m/sec.  turns  out  to  be 
4.83xl03  kg/cm2. 

From  the  results  which  were  obtained  for  wet  sand  and 
for  dry  sand,  we  see  that  the  wave  velocity  in  dry  sand  was 
twice  as  large  as  that  in  wet  sand,  and  Young's  modulus  in  dry 
sand  is  four  times  larger  than  for  wet  sand.  These  differences 
In  the  dynamical  behavior  between  dry  sand  and  wet  sand  are  in 
accordance  with  the  differences  observed  between  their  stiff¬ 
nesses  when  these  materials  are  handled  quasistatically . 
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5 .  Discussion  of  Results 

5.1.  General  Remarks 

In  the  previous  chater  It  was  mentioned,  that  the 
results  which  could  be  derived  from  the  experiments  by  using 
the  rigorous  analysis  presented  in  chapter  3,  were  very  limited. 
The  reason  for  this  is  that  small  changes  in  the  location  of 
the  coordinate  system  in  any  of  the  pulses,  as  well  as  small 
changes  in  the  shapes  of  the  pulses,  could  cause  enormous  errors 
in  the  computations  of  the  material  properties  of  the  lower  bar. 
This  is  proved  in  the  following  sections. 

5.2.  Effect  of  Coordinate  Shift  on  tan  6  *  and  C,-, 

Let  us  consider  an  incident  and  reflected  pulse  as 
shown  in  Fig.  33-  In  this  figure,  the  location  of  the  origin  of 
the  reflected  pulse  which  corresponds  to  the  origin  of  the 
Incident  pulse.  Is  not  welx  defined.  Hence  we  denote  by  0  the 
true  origin  of  the  reflected  pulse  (corresponding  to  the  e  ,  t 
coordinate  system),  and  by  O'  the  erratic  origin  of  this  pulse 
(corresponding  to  the  e  ,  t'  coordinate  system).  The  shift 
between  the  true  and  postulated  coordinate  systems  is  denoted 
by  m,  and  it  can  have  either  a  positive  or  a  negative  value . 

If  we  now  describe  the  Incident  pulse  as  a  Fourier  sum 
in  its  coordinate  system,  we  have: 

oo  oo 

e,  (t  =  a  +  l  a.  cos  kp  t  +  £  a.  sin  kp  t 

-L  °  k=l  KG  °  k=l  Ks  ° 


where : 
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a  =  —  f  e.(t)dt 
o  T  J  1 

o 


2 

=  m  e. ( t ) cos  kp  t  dt 
T  J  l  o 


=  ei(t)sin  kpQt  dt 

✓ 


Describing  now  the  reflected  pulse  in  the  postulated  coordinate 
system  (which  we  mistakingly  consider  to  be  the  true  coordinate 
system),  we  have: 

t  f  °°  T  °°  f 

e  (t1)  =  b  +  £  b,  cos  kp  t  +  £  b,  sin  kp  t 

r  o  k=1  *c  “o  k=1  ks 


where : 


-  J  er(t ' )  dt ' 


T 

2  f  ' 


ip  J  er(t  1  )cos  kpQt ’dt  ' 


T 

2  f  ’ 


ip  I  er(t  * ) sin  kpot ’ dt ' 


Describing  the  reflected  pulse  in  its  true  coordinate  system. 


we  have 


er(t)  =  bo  +  Z  bkc  COS  kpot  +  1  bks  sln  kpnb 


where : 


■  •  v&d*  *.n wkmnmmni l 
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m 

bo  =  T  er(t)dt 
o 

T 

2 

bkc  =  T  er(t)cos  kpQt  dt 

0 

m 

bks  =  I  er(t)sin  kpQt  dt 
o 

Since  the  reflected  pulse  is  practically  zero  in  the  small 
regions  near  its  edges,  we  can  also  write  its  true  Fourier 
coefficients  as  follows: 

T+m 

bo  =  f  f  er(t)dt 
m 

T+m 

2 

bkc  ~  T  I  £r(t)cos  kpQt  dt 
m 

T+m 

bks  -  |  er(t)sin  kpQt  dt 
m 

Now,  the  relation  between  uhe  mathematical  descriptions  of  the 
reflected  pulse  in  its  two  coordinate  systems  is: 

er(t)  =  er  ^ t_m) 

Hence,  the  true  Fourier  coefficients  of  the  reflected  pulse  can 


be  written  as: 
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T+m 

|  J  er(t-m)dt 

m 

T+m 

2  * 

^  er(t-m)cos  kpQt  dt 
m 

T+m 

p  f  • 

ip  j  er(t-m)sin  kpQt  dt 
m 


Substituting  the  relations  t-m  =  t'  and  dt  =  dt ’  in  the  last 
integrals,  and  changing  their  limits  accordingly,  wo  get: 


'ks 


o 


T 

If 


£r(t’)dt' 

er(t*)cos  kpo(t'+m)dt' 


er(t')[cos  kpQt ’  cos  kp^m-sin  kpQt ' 
er(t')sin  kpo(t*+m)dt' 

f 

er(t')[sin  kpQt '  cos  kpQm+cos  kpQt * 


o 


sin  kpQm]dt  ’ 


sin  kp  m]dt  ’ 
o 


Referring  now  back  to  the  expressions  for  the  Fourier 
coefficients  of  the  postulated  pulse  we  get  from  the  last 
expression  the  following  relations  between  the  true  and  the 
postulated  Fourier  coefficients: 


mm* 


I 


=  (cos  kpom)bkc  -  (sin  kpQm)bks 


(5D 


=  (cos  kpom)bks  +  (sin  kpom)bk( 


In  order  to  demonstrate  the  effect  of  the  shift  "m"  on  the 
computations  of  tan  6'  and  C2>  we  will  assume  for  simplicity  that 
the  upper  bar  is  perfectly  elastic,  that  is  tan  6=0.  In  this 
case  tan  6*  and  are  calculated  from  expressions  (21).  If  we 

f 

denote  by  (tan  6’)'  and  C2  the  erratic  properties  corresponding 
to  the  postulated  Fourier  coefficients,  and  by  tan  6'  and  the 
true  properties  corresponding  to  the  true  Fourier  coefficients, 
we  get  from  expressions  (21)  and  (51)  the  following  relations: 


(a)  ( tan  6 ' )  = 


-4(akcbks-aksb(c) 

a,2  +a2  -b!2-b;2 
kc  ks  kc.  ks 


.  ,pici/aLaakVbkc-bks>  _ _ 

°2  (akc+akS+bkc+bks)-2(ak-b)-+ak«b^> 


kc  kc  ks  kc 


(c)  tan  6 


•11(akcbks-aksbkc)coskpo,n-',<akcbkc+aksbks)sinkpo!" 

2  2  '2  '2 
a,  +a,  -o,  -b, 

kc  ks  kc  ks 


P  2^1 
Co  =  (r1-1) 


,2  2  .  '2  .  !2 , 
kc  aks  kc  ks 


(akc+aks+bkc+bks)-2(akcbkc+aksbks)coskpom+2(akcbks-aksbkc)sinkpc 
The  ratios  between  the  postulated  and  the  true  material 


properties  will  be: 


i 


t>5 


(tan  6 1 )  1 
tan  6  ' 


akcbks“aksbkc 


(akcbks-aksbkc)cos  kpom+<akc°kc+aksbks)sln  kpt 


m 


c. 


(53) 


,  P  p  I  0  *  p  t  t 

(a,  +a ,  +b,  ‘~+b,  )-2(a,  b,  +a,  b,  ) 

kc  ks  kc  ks'  kckc  ks  ks' 

If  we  assume  for  simplicity  that  the  Incident  pulse  is  a  symmetric 
pulse,  then  a^g  =  0,  and  relations  (53)  reduce  to: 


(a)  i^an  <'>'  » 

tan  6 


bkc 

cos  kp  m  +  —V—  sin  kp  m 
0  '  o 

bks 


(b)  ^  - 


(5*0 


(ali+bk^^-2(akcbk<:)c°s  kpom  +  2(akc'ks)5ln  kpom 

<*L«>  -  ^kcVf 


Now,  since  the  .incident  pulse  was  assumed  to  be  symmetric,  the 
ieflected  pulse  will  be  very  close  to  a  symmetric  pulse.  The 
reason  for  this  is,  as  will  be  shown  later,  that  for  small  tan  6' 
the  reflected  pulse  has  almost  the  same  shape  as  the  incident 

>2  p  i  o 

pulse.  Hence,  b.  can  be  neglected  in  comparison  to  a,  and  b.  , 

kc  xc  ’ 

and  expression  (5^b)  becomes  after  some  algebraic  manipulations: 


,  kp  m 

o 


C2  _  ,  ,  „  akc(bkc  tan  — 

('  Id  to 

2  (akc'bkc) 


+  bks> 


sin  kp  in 
o 


(5^c  ) 
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( tan  <5  * ) 1 

From  expression  (5^a)  we  see  that  Arr— fi  is  strongly  depend- 

u  an  o 

ent  on  the  shift  "m"  even  for  small  shifts.  The  reason  is  that 

4 


bkc 

the  coefficient  r~r^  is  very  large  (because  bVc,  is  very  small  in 


ks 

f 

comparison  to  bkc,  as  was  mentioned  above),  and  hence  the  second 

term  in  the  denominator  of  expression  (5^a)  may  become  equal  or 

larger  (in  absolute  value)  than  the  first  term  in  the  denominator 

of  that  expression.  , 

bkc 

If  we  assume  that  — —  is  positive,  then  for  positive  shifts 

bks 

( 6*)' 

m,  ~fra"— g  4 —  may  become  much  smaller  than  1.0.  For  negative  m, 

( fc 3.n  6  1  )  ’ 

tarT  6  r~  may  136001116  much  larger  than  1.0  which  may  cause  very 
large  errors  in  the  value  of  tan  6’.  In  fact  when  m  is  negative, 
the  whole  denominator  in  (5^a)  may  turn  out  to  be  negative,  which 
leads  to  a  negative  (tan  6')’.  Hence  we  see  that  relatively  small 
shifts  (of  the  order  of  1%  of  the  pulse  width)  may  lead  to  absurd 
results  for  tan  6  * .  , 

c2 

From  expression  (5^c)  we  see  that  is  moderately  depend- 

b2 

ent  on  m.  The  reason  is  that  the  second  term  in  the  right  hand 
side  of  expression  (5^c)  is  much  smaller  than  1.0  at  least  for 


the  first  few  harmonics 


r  * 

b  2 


Hence  is  very  close  to  1.0  for  the 

L  2  kp  m 

first  few  harmonics.  However  for  the  higher  harmonics,  tan  -0-°  - 

C  '  d 

2 

may  become  very  large,  causing  large  errors  in  .  As  will  be 

b2 

shewn  later,  the  effect  of  small  changes  in  the  pulse  shape  may 
introduce  large  errors  in  C.;  even  at  the  lower  harmonics  (with 
the  exception  of  the  first).  Hence  the  only  value  of  which 
can  be  obtained  with  reasonable  accuracy  is  uhe  value  of  C0  for 


am.'- ' 


the  first  harmonic.  All  the  arguments  which  were  presented  till 
now,  are  still  valid  when  the  incident  pulse  has  an  arbitrary 
shape  and  not  necessarily  a  symmetric  one.  In  this  case  all  the 
Fourier  coefficients  in  expressions  (53)  must  be  taken  into 
account,  and  the  analysis  becomes  somewhat  more  involved.  How¬ 
ever,  the  conclusions  which  may  be  derived  in  this  case  are  the 
same  as  those  which  were  obtained  for  a  symmetric  incident  pulse. 
In  fact,  when  the  pulses  of  Fig.  12b  (which  are  all  asymmetric) 
were  used  in  the  computations,  the  resulting  values  for  tan  6' 
at  the  first  harmonic  varied  between  0.6  and  -0.5  for  shifts 
between  +2%  and  -2%  of  the  pulse  width  (the  true  value  for  tan  6' 
in  this  case  is  0.1).  Similar  cr  larger  variations  in  tan  6' 
were  also  obtained  for  the  higher  harmonics.  The  resulting 
values  for  C 2  which  were  obtained  from  Fig.  12b  varied  between 
890  m/sec  and  910  m/sec.  for  the  first  harmonic  (about  2%  error) 
due  to  the  shift  effect. 

The  values  of  C2  for  the  higher  harmonics  differed  con¬ 
siderably  from  these  values  because  of  the  shape  effect  which 
will  be  discussed  later.  However  the  deviations  of  C2  at  each 
harmonic  from  the  mean  value  of  C2  at  this  harmonic,  increased 
gradually  from  +1%  at  the  first  harmonic  to  +50%  at  the  tenth 
harmonic  due  to  shifts  of  +_2%  of  the  pulse  width. 

5 . 3  Effect  of  Coordinate  Shift  on  E2(t) 

Suppose  that  the  incident  and  the  reflected  pulses  are  again 
those  shown  In  Fig.  33.  Now  the  stress  relaxation  function  E0(t) 
may  be  computed  from  expressions  (47)  and  (49)  which  yield 
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E2(t)  -  [pE2]  _Q<5  = 


(plCl)2r/ei+er/. 


e .  — £ 
i  r 


From  a  physical  point  of  view,  it  is  to  be  expected  that  E2(t) 
should  be  a  slowly  decreasing  function  w.r.t.  time.  Now,  the 
time  t  is  inversely  proportional  to  the  transform  variable 
p(p  =  .  Thus  E2(t)  (or  pE2(p))  must  be  a  slowly  increasing 

function  w.r.t.  p.  We  therefore  conclude  that  the  slope 
d(pEp)  dE0(t) 

— - — —  (or  — 5 - )  must  be  positive  and  close  to  zero. 

dp  '  dp 

Suppose  now  that  the  origin  of  the  reflected  pulse  was 
shifted  a  distance  "m"  from  its  exact  location.  The  Laplace 
transform  of  the  reflected  pulse  w.r.t.  to  its  erratic  coordinate 
system  will  be  e-™13^.  Inserting  this  value  of  the  transformed 

I 

reflected  pulse  in  expression  (55) ,  and  denoting  by  E2(t)  the 
resulting  erroneous  stress  relaxation  function,  we  get: 


E2(t) 


(p  C  )2  e±+e  pe 

- — — — [(-i - -)  ] 

P0  T  »-mP7  OR 

^  Ei_e  br  p=^-2. 

Ty 


(56) 


Suppose  now  that  the  shift  m  is  very  small.  In  this  case  we  can 
write  approximately: 


Substituting  this  value  for  e~mp  in  expression  (56),  performing 
some  algebraic  manipulations  and  retaining  only  first  order  terms 


in  m,  we  get: 


^P1°V~  r/ei+cr^  ’  „  ei£r  ^  ei+er  ^  n 

Ep  ( t )  -  —  [  (z — =~)  -  **mp  —z— 3 — 3 3 


£i"£r 


(ei-£r)' 


The  erroneous  slope 


dE2(t) 


-  will  be : 


-  ta  uP 

dp  p2  ldp  dp  (7^)3 


The  first  term  in  the  right  hand  side  of  the  last  expression  is 
simply  the  slope  of  the  exact  stress  relaxation  function  Ep(t). 
Thus  we  have: 

dEj(t)  dE2(t)  (plCl)2  d 

■TT  *  —  -  «-  — —  5F[P  (57) 

dE2(t) 

Now,  since  — ^ —  is  small,  the  second  term  in  the  right  hand 
side  of  expression  (57)  may  be  equal  or  larger  (in  absolute 
value)  than  the  first  term  in  that  expression,  even  when  "m;i  is 


very  small.  The  second  term  may  also  have  either  positive  or 


negative  values  depending  on  the  sign  of  m,  and  the  sign  of 

_  __  _  f 

j  e  j  £p  ( £ )  dE p  ( t ) 

— r — r p - 5 — ].  Hence  — -s -  may  come  out  as  many  times 

dp  (7  7  \3  dp 

v  i  r; 


larger  or  many  times  smaller  than 


dEp ( t ) 


,  cr  may  even  assume  a 


negative  value  . 


In  fact,  when  Ep(t)  was  computed  from  Pig.  12b,  a  slowly 
increasing  function  w.r.t.  p  (or  slowly  decreasing  function  w.r.t 


t)  was  obtained  for  zero  shift.  However  the  slope  of  this  func¬ 
tion  was  tripled  when  a  shift  of  +25?  was  introduced,  and  turned 


negative  when  a  shift  of  -2%  was  Introduced  (a  ne gar  cpe 

means  a  decreasing  function  w.r.t.  p  or  an  increase  tlon 

w. r. t.  t ) . 

Thus  we  see  that  small  errors  in  the  location  of  the 
coordinate  system  In  the  reflected  pulse,  may  introduce  large 
errors  in  the  slope  of  the  stress  relaxation  function  E,-,(t),  and 
hence  cause  the  curve  of  E2(t)  vs,  t  to  slope  in  the  wrong 
direction. 


5 . 4  Effect  of  Pulse  Shape  on  tan  6  *  and  C 2 

In  order  to  investigate  the  effect  of  small  changes  in  the 
pulse  shapes  on  tan  6',  let  us  assume  that  two  similar  experi¬ 
ments  gave  exactly  the  same  incident  pulses  but  slightly  differ¬ 
ent  reflected  pulses.  We  will  again  assume  for  simplicity  that 


the  incident  pulses  are  symmetric  (yielding  a^g  f  0)  and  hence 
the  reflected  pulses  are  very  close  to  being  symmetric.  Denoting 
by  the  subscript  (1)  all  the  quantities  which  belong  to  the  first 
reflected  pulse,  and  by  the  subscript  (2)  all  the  quantities 


which  belong  to  the  second  reflected  pulse,  we  have  from  the 
first  expression  In  (21): 


tan  6 


^akcbks  ( 1 ) 

2  2  .2 
ake-Dkc ( 1 )  ’  ks ( 1 ) 


tan  <5 


(2) 


^akcbks(2) 

2  2  2 
akc-bkc(2)~bks(2) 


t  T 

Dividing  tan  <5  ( -jj  by  tan  6^^  we  get: 
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tan  S(’l)  _  ,ako-bkc(2)-bks(2) 
»  “  v  2 
tan  <5(2) 


_ 

ate-bk=(i)--bL(i)  bks<2) 


(58) 


Now  and  b^s(2)  are  bobh  very  small  quantities  in  compari¬ 

son  to  all  the  other  Fourier  coefficients  (for  tan  6 ’  =  0  we 
will  get  an  exactly  symmetric  reflected  pulse,  and  hence  b^s  =  0). 
However  the  ratio  between  them  may  be  much  larger  than  unity  or 
much  smaller  than  unity  or  even  negative. 

Hence  we  see  that  small  changes  in  the  pulse  shape  may 
lead  to  completely  different  results  for  tan  6*.  In  order  to 
illustrate  all  the  above  considerations,  let  us  assume  that  the 
incident  pulse  has  the  form  of  a  half  sine  wave  in  the  interval 
0 — 7t ,  and  the  reflected  pulses  will  then  have  the  form  of  a 
slightly  asymmetric  half  sine  wave  in  the  same  interval. 

The  suitable  expressions  for  all  these  pulses  may  be 
written  as : 


e.  =  A  sin  t 

l 

a,  t 

=  Be  sin  t 

a?t 

er(2)  =  Be  sin  t 

Here  and  are  small  numbers  and  determine  the  .location  of 
tne  pulse  peak  (which  is  very  close  to  tt/2). 

Computing  the  Fourier  coefficients  of  all  three  pulses 

2 

and  neglecting  a  w.r.t.  unity  we  get: 
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bkc(l) 
bkc ( 2 ) 

bks(i) 

bks(2) 


4  A 


r r  4k  -1 

a,  it 

2 { e  1  +1)  B 
11  4k2-l 

a  tt 

2 (e  ^  41)  B 
11  4k2-l 


a1  tt 


8a  (e  "  +1)  kB 


TT 


(4k2-!)2 


a  tt 

8a^(e  +1) 


kB 


*  (4k2-!)2 


Since  and  a2  are  small  numbers  (say  G.04  and  0.02)  we  set: 


bkc(l)  ~  bkc ( 2 ) 


4B  1 


4k  -1 


l6a 


ks ( 1 ) ' 


ks  ( 2 ) 


1  kB 


"  (Hk2-1)2 

l6a. 


2  kB 


"  (Kk2-!)2 


o  2 

Thus,  from  expression  (?8)  and  by  neglecting  and 

w.r.t.  all  the  other  Fourier  coefficients,  we  get: 

tan  6(1)  _  ^1 
tan  <5^2)  ^ 


Now  an  can  be  two  or  three  times  larger  than  or  even  have  the 

I 

opposite  sign  of  a?.  Hence  tan  5^^  may  differ  considerably  from 

I 

tan  <5(2)  anc*  this  results  from  Just  small  changes  in  the  shape  of 
the  reflected  pulse. 


In  order  to  investigate  the  effect  of  small  changes  in 
the  pulse  shapes  on  C^j  let  us  assume  now  that  two  similar  experi¬ 
ments  result  in  the  same  reflected  pulses  but  from  two  slightly 
different  incident  pulses  which  had  the  same  maximum  amplitudes. 

We  will  again  assume  for  simplicity  that  the  incident  pulses  are 
symmetric.  Denoting  by  the  subscript  (1)  all  the  quantities 
belonging  to  the  first  incident  pulse*  and  by  the  subscript  (2) 
all  the  quantities  which  belong  to  the  second  Incident  pulse,  we 
get  from  the  second  expression  in  (21): 


'2(1) 


Pl^l 

--  (-^) 
P2 


2  2  2 
akc  ( 1 )  ~bkc  ~bks 

(akc(l)'bkc)  +bks 


C2(2) 


P1C1 

p2 


2  2  2 
akc  ( 2 )  ~bkc~bks 

(akc(2)"bkc)  +bks 


2 

Dividing  C2q)  by  c2(2)  and  neglecting  bks  in  comparison  to  all 
the  other  Fourier  coefficients,  we  get: 


C  2  ( I ) 
C  2  ( 2  ) 


/akc  (1  )~l~bkc  s 

akc  ( 1 )  ~bkc  =  |akc(l)+bkc^akc(2)~bkc^ 
( akc ( 2 )4bkc s  ^  akc ( 1 ) -bkc  ^  ^  akc ( 2 )+bkc^ 
akc ( 2 )  ~bkc 


(59) 


Now  suppose  that  the  first  incident  pulse  is  a  half  sine  wave  in 
the  interval  0  —  tt ,  and  the  second  incident  pulse  is  a  parabola  in 
the  same  Interval.  Those  two  pulses,  when  plotted,  are  very 
similar  to  each  other.  The  reflected  pulse  will  be  assumed  as 
before  to  have  the  shape  of  a  slightly  asymmetric  half  -ine  wave. 
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The  suitable  expressions  for  all  three  pulses  are: 

ei(l)  =  A  sln 

_  11 A  ,  .  ,2s 

£i(2)  ~~  ~i 2  ' 

er  =  Be0^  sin  t 

The  Fourier  coefficients  of  the  three  pulses  are: 


'kc(l) 

4A  1 

’  77  4k2-l 

kc  (2 ) 

4A  1 

"  TT2  k2 

ks  (1 ) 

aks (2)  =  0 

2(eaU+l)  B 

4B  .1 

(for  small  a) 

kc 

*  Hk2-1 

17  4k.2-l 

8ct(eaTr+l)  Bk 

l6aB 

k 

ks 

77  (4k2-!)2 

TT 

(4k2-!)2 

(for  small  a) . 


Substituting  the  Fourier  coefficients  in  expression  (59)  we  get 


'2(1) 

'2(2) 


4  A 


4k  -l 


4b 


TT 


4k2-l 


)( 


4A  1_ 

TT2  k2 


4B 

TT 


4  k  2  - 1 


4  A 


4k  -1 


4B  1  w  4A  1  4b 

■„  2  ,  M  2  .  2 
TT  4  k  - 1  TT  k  TT 


4k.2-l 


or 


C2(l) 
C  2  (  2  ) 


(1  +  |)( 


i 


+ 


B 

A 

¥ 

A 


4k2-l 
_ 1 

4k2-l 


) 


) 
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Now  let  us  suppose  for  example  that  the  amplitude  ratio  B/A  is 
equal  to  0.5.  Computing  c2(i)/,C2(2)  for  the  flrst  three  harmonics 
(k=l,  k-2,  k=3)  we  get: 


k=l 

k=2 

k=3 

0.935 

1.230 

1.270 

[2(1) 

'liO. 


Thus  we  see  that  small  changes  in  the  shape  of  the  incident  pulse, 
though  not  effecting  the  value  of  C 2  at  the  first  harmonic  too 
much  ( 7 %  error  in  our  example),  may  introduce  large  errors  in 
the  values  of  at  the  higher  harmonics  (in  our  example  we  have 

23#  error  in  C2  of  the  second  harmonic  and  27#  error  in  C 2  of 
the  third  harmonic).  The  reason  for  the  sensitivity  of  the  value 
of  C 2  (of  the  higher  harmonics)  to  small  changes  in  the  pulse 
shape  is,  that  the  Fourier  coeff icients ,  excluding  those  of  the 
first  harmonic,  are  very  dependent  on  the  exact  shape  of  the 
pulse.  Thus,  two  pulses  which  have  almost  similar  shapes,  may 
have  entirely  different  higher  Fourier  coefficients.  The 
conclusions  of  the  discussion  presented  up  to  now  are  quite 
general  and  will  hold  for  pulses  of  arbitrary  shape  when  subject 
to  small  changes  in  their  shapes.  When  the  pulses  of  Fig.  12b 
were  analyzed,  the  resulting  values  for  tan  5 ’  for  the  first  ten 
harmonics  varied  between  -in  i  and  1 . k  for  a  specific  location 
of  the  coordinate  syste  ^.n  the  reflected  pulse. 

(This  location  was  the  one  which  gave  the  correct  result 
tan  6'  =  0.1  at  the  first  harmonic.  Other  locations  of  the 
coordinate  system  resulted  in  even  larger  variations  in  tan  6  ' 
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computed  for  different  harmonics.)  The  large  variations  in 
tan  6'  are  attributed  mainly  to  small  Inaccuracies  in  the  shapes 
of  the  pulses  in  Pig.  12b.  These  inaccuracies  are  probably 
Introduced  by  small  experimental  errors.  When  the  phase  velocity 
Cg  was  computed  from  Fig.  12b  for  the  first  ten  harmonics,  values 
between  600  m/sec.  and  1770  m/sec.  were  obtained.  (These  values 
were  obtained  for  the  same  location  of  the  coordinate  system 
which  was  mentioned  oefore.  For  other  locations  of  the  coordin¬ 
ate  system  larger  variations  in  were  obtained.)  The  value 
which  was  obtained  for  C2  for  the  first  harmonic  was  900  m/sec. 
which  is  probably  close  to  the  correct  value  for  this  harmonic. 

(A  similar  vaxue  for  C-,  for  the  first  harmonic  was  also  obtained 

£ 

from  Fig .  12a. ) 

From  the  discussion  presented  in  this  paragraph,  we  see 
that  the  computed  values  for  tan  6 '  and  C 2  (for  the  higher 
harmonics)  are  very  sensitive  to  small  changes  in  the  pulse 
shape.  The  only  quantity  which  is  not  effected  too  much  by 
these  small  changes  and  hence  could  be  computed  with  reasonable 
accuracy,  is  the  phase  velocity  C2  for  the  first  harmonic. 

5 ..  5  Dependence  of  the  Pulse  Shape  on  tan  6  ' 

In  the  previous  paragraphs  it  was  mentioned  that  for  small 
tan  o',  the  shape  of  the  reflected  pulse  is  almost  similar  to 
the  shape  of  the  incident  pulse  (although  the  amplitudes  are  in 
general  different.)  A  proof  of  this  statement  will  now  be  given. 


a..iX 
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Before  starting  on  this  proof,  however,  we  need  to  know 
the  relations  between  the  Fourier  coefficients  of  the  incident 
strain  pulse  and  those  of  the  reflected,  strain  pulse.  This  can 
be  done  by  substituting  relations  (19)  in  the  first  two  relations 
in  (8)  yielding  (after  some  algebraic  manipulations): 


,  e  2s ( ) 

_  (  -LzS  ■)  a  4- _ ^  _  _  o 

1+8  kc  (l+s)2  ks 


ll+s;  ks 


2s(62-61) 

(l+s)^ 


(60) 


Suppose  now  that  the  incident  pulse  is  given  in  the  basic  interval 


0-2tt  by : 


e±  =  0.0687eO,19(2,4“tM';?(27r-t)1,2 


This  pulse  is  shown  in  Fig.  34.  Suppose  also  that  C 2  is  changing 
with  frequency  as  predicted  by  the  linear  theory  of  viscoelasticity 
for  constant  tan  6  which  yields: 


C2  ~  C20('1  + 


tar^Ll  £n 
75 


(61) 


Here  C'2Q  is  the  phase  velocity  at  the  basic  frequency  p  .  Our 
object  now  is  to  predict  the  shape  of  the  reflected  pulse  for 
different  values  of  tan  6’,  assuming  that  C2  is  varying  with 
frequency  according  to  expression  (61).  We  will  also  assume  for 
simplicity  that  the  first  bar  is  elastic.  The  first  step  is  to 
find  the  Fourier  coefficients  of  the  incident  pulse.  From  these 
coefficients,  the  corresponding  Fourier  coefficients  of  the 
reflected  pulse  may  be  found  uy  using  expressions  (60)  and  (61). 
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(Expression  (6l)  is  used  to  compute  "S"  for  each  frequency.) 

From  the  Fourier  coefficients  of  the  reflected  pulse,  the  shape 
of  this  pulse  can  be  constructed. 

The  procedure  just  described  was  carried  on  with  the  aid 
of  the  computer,  for  two  values  of  ** S'*  at  the  basic  frequency 
(0.5  and  0.75)  and  for  five  values  of  tan  6'  (0,  0.05,  0.1,  0.2, 
0.3).  The  results  are  shown  in  Figs.  35a,  35b,  35c,  35d. 

Figures  35a  and  35c  give  the  immediate  results  of  the  computa¬ 
tions,  from  which  it  might  appear  (at  least  from  Fig.  35c)  that 
there  are  considerable  differences  between  pulse  shapes  corre¬ 
sponding  to  different  values  of  tan  6'.  However  a  closer  look 
at  these  figures  reveals  that  in  fact  all  the  pulses  are  almost 
identical,  and  they  are  just  shifted  a  small  distance  from  each 
other.  Thus,  if  all  the  pulses  are  shifted  back  towards  each 
other,  the  results  which  are  obtained  are  those  shown  in  Figs  . 

35b  and  35d.  From  Figs.  35b  and  35d  we  see  that  any  small  value 
which  is  assumed  for  tan  6 T  results  in  almost  the  same  reflected 
pulse.  Thus  it  is  again  verified  that  tan  6'  cannot  be  derived 
with  sufficient  accuracy  from  the  strain  pulses  which  v ere 
obtained  in  the  experiments  described  in  this  report. 

It  is  interesting  to  see  from  Figs .  35b  and  35d  that  the 
reflected  pulses  for  different  values  of  tan  6',  almost  coincide 
with  the  reflected  pulse  for  tan  6 ’  =  0  which  corresponds  to  a 
perfectly  elastic  second  bar.  This  provides  some  indication  that 
the  phase  velocity  for  the  first  ha  .nonic  can  be  approximately 
computed  by  assuming  tnat  the  two  >.  irs  are  elastic,  and  by 
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measuring  the  maximum  amplitudes  of  the  pulses  and  inserting 
them  in  expression  (50)  which  in  turn  yields  a  straightforward 
result  for  Cp • 

From  expressions  (5)  which  correspond  to  the  case  of  two 
elastic  bars  attached  together  we  have: 


pC)C„~p1  Ch 

Er<‘>  '  (p&p  e-  > 


Thus  we  see  that  in  this  case,  the  shapes  of  the  incident  and 
the  reflected  pulses  are  similar  (although  the  amplitudes  of  the 
pulses  are  different).  Now,  since  the  shapes  of  the  reflected 
pulses  when  tan  6 »  /  0  are  very  close  to  the  shape  of  the 
reflected  pulse  when  tan  S’  =  0,  we  conclude  that  the  reflected 
pulses  are  always  close  in  shape  to  the  incident  pulse  when 
tan  6'  is  sufficiently  small  (say  less  than  0.3,  which  includes 
all  the  conventional  plastics).  The  conclusions  present  in  this 
paragraph  are  quite  general  and  were  also  derived  when  different 
incident  pulses  were  assumed. 


5 . 6  Effect  of  Pulse  Shape  on  E ^ ( t ) 

In  order  to  investigate  the  effect  of  small  changes  in  the 
pulse  shape  on  Ep(t)  we  will  assume  that  Ep(t)  may  be  approxi¬ 
mated  to  in  a  limited  interval  by: 

Ep ( t )  =  Eoe_at  (62) 

For  the  time  period  0-200  poec  (which  Is  the  approximate  duration 
of  most  of  the  strain  pulses  obtained  in  the  experiments),  one 
may  expect  a  reduction  of  about  c>%  in  K0(t)  .  For  illustrative 
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purposes  however,  we  will  assume  changes  in  E2(t)  which  will  vary 
between  10%  reduction  and  5%  increase.  Thus  we  may  write: 

E2(2Q0)  =  Eoe_Ct,20°  =  3Eq 

where  0.90  <_  3  <_  1.05. 

Prom  the  last  expression  we  get: 


Hence : 


a 


Jin  3 
200 


8  =  1.05  -  a  = 


-  2.439*10“^  y sec-1 


8  =  1.0  a  =  0 

8  =  0.9  a  =  +5. 26  8x10“^  ysec 


The  graphs  E2(t)  vs.  t  fcr  these  values  of  a  are  shown  in  Fig.  36. 
Now,  assuming  for  E2(t)  the  form  given  in  expression  (62),  we 
may  use  expression  (46)  in  order  to  determine  the  relation  between 
the  shapes  of  the  incident  and  the  reflected  strain  pulses.  Doing 
this  we  get : 


and  retain  only  first  order  terms  in  a  .  Here  we  also  assume 

that  p  is  large  compared  with  a  so  that  -  is  still  much  less 

than  unity.  If  p  becomes  small,  which  corresponds  to  t  becoming 

0  5 

large  (recall  p  -  —^L) ,  the  Inaccuracy  In  the  mathematical 
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1 

procedure  gradually  Increases.  Evaluating  — —  in  a  Taylor 

1+- 

P 

series  and  performing  some  other  algebraic  manipulations  in  (63), 
in  which  we  only  retain  first  order  terms  in  a,  we  finally  get: 


-  _r^2Eo  "  P1C1 

a  ={x- 


plClf2Ec 


a  1  — 

2  ^  Gi 


fK  +  +  PjCp  P 

Taking  the  inverse  Laplace  transform  of  the  last  expression  we 
get : 

t 


/4.\  _  /p2Eo  "  P1C1  P1C1  /p2Ec 

e  (t)  =  *- — — -  e,(t) 

/  n  .  _  r\  -L 


v/fTEo  +  fl0! 


(/P2Eo  +  Pi0!5 


a 

< 


(r)dT  (64) 


Expression  (64)  gives  us  the  approximate  relation  between  er(t) 
and  e^Ct)  for  times  which  are  not  too  large  (say  not  larger  than 
the  pulse  duration  which  is  200  psec).  In  order  to  Illustrate 
the  dependence  of  er(t)  on  a,  let  us  assume  that  e^(t)  Is  given 
as : 

0  <  t  <  T0 
*  >  To 

Here  Tq  is  the  duration  of  the  incident  pulse  which  Is 
assumed  to  be  200  psec.  Computing  er(t)  from  expression  (64)  we 
get : 


/p2Eo 

Assuming  for  example  that  -—-77-  ■  =  0.5*s  and  plotting  e  (t)  for 

P1C1  r  , 

the  three  values  of  a  which  were  given  above  (ot=-2 . 439*10  ,  a=0, 
a=+5. 268x10  S  we  get  the  pulse  shapes  which  are  shown  in  Fig.  37a. 
From  Fig.  37a  it  seems  that  there  are  some  differences  between 


the  pulses  er(t)  which  correspond  to  different  values  of  a  , 
especially  near  the  tails  of  the  pulses. 


However,  the  mathematical  treatment  which  results  in  Fig. 37a 
is  not  accurate  for  large  t,  and  in  fact  all  the  reflected  pulses 
almost  coincide  also  at  the  ends  of  their  periods.  This  may  be 
proved  with  the  aid  of  a  theorem,  which  relates  the  behavior  of 
a  function  of  t  at  t  +  ®  to  the  behavior  of  its  Laplace  transform 
when  p  -*■  0 .  This  theorem  states: 


11m  e  (t)  =  lim  p  e 
t->°°  p+o 


If  we  define  a  characteristic  pulse  velocity  C?  in  the  second  bar, 

2  2Eo  PpCp 

such  that  E  -  PpCo  then,  — —  p  =  —  p —  which  is  the  ratio  of  some 

pi°i  piui 

typical  characteristic  impedance  in  the  second  bar  to  the  charac¬ 
teristic  impedance  of  the  first  bar.  A  typical  value  for  this 
ratio  may  be  taken  as  0.6. 
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Applying  this  theorem  to  expression  (63)  (which  is  valid  for  any 
value  of  a,p,  and  t)  we  get: 


1 im  e  ( t ) 
t-*00  r> 


lim  p  er 
p->-o 


lim  p 
p+o 


-J- 

7 


V 

2p+a 


- 


E  p 
2  p+a 


Dici 


0  (65) 


The  last  is  true  because  e.  approaches  a  constant  value  when 

TT  t 

p  -*  0.  (In  fact  for  =  sin  —  at  t  £  Tq,  we  get 


£ .  = 


-pT 

(f-)(l+e  ° 

t~ j? 

o 


) 


from  which  lim  e^  =■ 
p+o 


(£-) 

o 


=  const . )  From 


(65) 


we  see  that  all  the  reflected  pulses  for  all  values  of  a  must 
approach  zero  when  t  is  sufficiently  large.  Taking  this  last 
fact  into  account,  it  seems  reasonable  that  all  the  reflected 
pulses  will  look  like  the  pulses  which  are  shown  in  Fig.  37b 
rather  than  those  which  are  shown  in  Fig.  37a. 

From  Fig.  37b  we  see  that  for  values  of  a  between 
-2.439x10-^  ysec  ^  and  +5.268x10”^  usee  \  which  corresponds 
to  changes  in  Eo(t)  between  a  5?  increase  and  a  10?  reduction, 
we  get  essentially  the  same  reflected  pulse.  Thus  we  see  that 
the  graph  of  E0(t)  versus  t,  cannot  be  derived  with  sufficient 
accuracy  from  the  experimental  results.  In  other  words,  very 
small  changes  in  the  shape  of  the  reflected  pulse  may  lead  to 
large  errors  in  the  determination  of  the  trend  of  the  graph 
Ep(t)  versus  t.  However  the  nve-age  value  of  En(t)  in  the 
limited  interval  of  time  0-200  ysec,  may  be  taken  as  the  constant 
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value  Eq  corresponding  to  a=Q.  This  value  was  termed  the  effec¬ 
tive  dynamic  Young's  modulus  in  the  previous  chapter,  and  since 
it  corresponds  to  a  perfectly  elastic  second  bar  (a=0  means 
E2(t)  =  const,  which  corresponds  to  an  elastic  material),  it  can 

p 

be  computed  from  the  relation  E2(t)  =  pC2  where  is  the  pulse 
velocity  in  the  second  bar. 


5 . 7  Concluding  Remarks 

It  has  been  shown  that  most  of  the  material  properties  of 
the  second  bar  cannot  be  derived  by  the  experimental  technique 
described  in  this  report. 

This  applies  especially  to  tan  5'  which  should  be  derivable 
from  the  first  expression  in  (21)  (assuming  that  the  first  bar 
is  perfectly  elastic).  The  value  which  is  computed  from  this 
expression  is  mainly  governed  by  the  numerator  a.^c s  — £3.^ sb^ ^  . 

If  both  bars  were  elastic,  the  incident  and  the  reflected  pulses 
/ould  have  had  exactly  the  same  shapes  so  that  e^(t)  =  Ae^(t) 
where  A  is  some  constant. 

In  this  case: 


b, 

=  A  a, 

kc 

kc 

b, 

=  A  3, 

ks 

ks 

and  hence : 


akc^ks  aki 


kc 


‘  akCUt*k! 


)-aksUakc)  " 


0. 


Thus  we  see  that  in  this  case  tan  6'  turns  out  to  be  aero  which 


was  expected  a  priori. 
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However,  when  tan  S'  /  0,  the  reflected  pulse  has  a  shape 

just  slightly  different  from  the  shape  of  the  incident  pulse. 

In  this  case  the  difference  a,  b,  -a,  b,  is  a  very  small 

kc  ks  ks  kc  J 

number,  and  hence  relatively  small  changes  in  any  of  the  Fourier 
coefficients  may  considerably  effect  this  difference,  which  in 
turn  might  result  in  enormous  errors  in  the  results  for  tan  6 ' . 
Now,  changes  in  the  Fourier  coefficients  may  be  caused  by  small 
errors  in  the  zero  of  the  coordinate  system  in  the  reflected 
pulse  (see  expressions  (51))  or  by  small  changes  in  the  shapes 
of  the  pulses.  Both  effects  are  always  present  in  any  experi¬ 
ment  because  of  inherent  experimental  errors. 

The  effects  of  small  changes  in  the  coordinate  system 
location  and  small  changes  in  the  pulse  shapes  on  C 2  and  E2(t), 
are  not  so  obvious  an  first  sight  as  they  are  in  the  case  of 
tan  6*.  However  some  explanation  and  some  examples  were  brought, 
which  showed  the  large  effect  of  small  inaccuracies  in  the 
coordinate  location  and  pulse  shapes  on  and  E0(t). 

In  the  discussion  given  in  this  chapter,  it  was  shown 
that  the  only  quantity  which  may  be  determined  from  the  experi¬ 
ments  with  reasonable  accuracy  is  the  phase  velocity  for  the 
basic  frequency  p  . 
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6 .  Summary 

An  experimental  technique  was  described  in  this  report  from 
which  theoretically,  the  material  properties  of  viscoelastic 
materials  could  be  derived.  However  from  the  experimental 
results  and  the  discussion  which  followed,  it  the 

only  quantity  which  could  be  determined  with  sufficient  accuracy 
was  the  phase  velocity  at  the  basic  frequency  in  the  viscoelastic 
material.  It  was  also  shown  that  in  order  to  determine  approxi¬ 
mately  this  phase  velocity,  there  is  no  need  to  use  a  rigorous 
analysis.  Instead,  direct  measurements  and  straightforward 
calculations  based  on  elastic  theory  give  the  desired  result  with 
less  effort  ar.a  the  same  degree  of  accuracy. 


87 


Appendix  I  -  Wave  Propagation  In  Viscoelastic  Bars 

Consider  an  extensional  wave  traveling  In  the  x  direction 
along  a  thin  viscoelastic  bar.  The  equation  of  motion  is: 


P 


da 

dx 


where  u  is  the  displacement  in  the  x  direction,  p  is  the  density 
of  the  bar  and  a  is  the  stress. 

In  order  to  solve  this  equation,  we  may  try  a  solution  of 
t  he  f o  rm : 

ip (t-  3) 

u  =  u  e  C 

o 


This  solution  corresponds  to  a  sinusoidal  wave  of  angular 
frequency  p  traveling  with  phase  velocity  C  along  the  bar  in  the 
direction  of  increasing  x.  From  the  expression  for  the  displace¬ 
ment  we  get : 


a 


<E1+1E2>  If 


-  (Ei+iE2)uoe 


ip  ( t- 


C 


where  E^+iE2  is  the  complex  Young's  modulus  of  the  material. 
Substituting  the  expressions  for  u  and  a  in  the  equation  of 
motion,  we  get: 

2  ip  i  t-  ~)  2  ip ( t - 

-  p  ^  u  e  C  =  -2_(E1+iE„)ue  ^ 

C2  0  C2  1  2  0 

from  which1* 


IT  _ 

We  see  that  C  turns  out  to  be  a  complex  phase  velocity. 


t 

t 
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«  E, +iE^ 
p-^  -  l  c-  _  _ 

P  P 


Ee 


16 


Here : 


-P 


1+E2 


tan  6  = 


Substituting  the  result  for  G  in  the  expression  for  u,  we  get: 


ip  ( t- 


u  =  u  e 
o 


/i 


E/p  e 


16/2 


=  uQe 


ip[t-  x(cos  6/2-i  sin  6/2 )  j 
V'e/p 


The  last  expression  may  be  written  as: 

ip ( t-  £) 


u  =  u  e”ax  e 
o 


where  C 


:  =  /e/ 


p  sec  —  is  the  phase  velocity  of  the  wave  and 


p  tan  j 

a  = - ^ - is  the  attentuation  coefficient. 

The  above  analysis  also  holds  for  a  wave  travel  lug  in  the 
direction  of  decreasing  x,  in  which  case  we  get: 


,  ip(t+  5/ 
tax  1  C 
u  =  u  e  e 
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dP 

Appendix  IV  -  Calculation  of  A™  for  Viscoelastic  Materials 

The  phase  velocity  in  a  linear  viscoelastic  material  with 
a  small  tan  6,  is  dependent  on  the  frequency  according  to  the 
following  relation  (see  Kolsky  [11]): 

(a)  C  -  C  (1  +  £n  E_) 

Here  C  is  the  phase  velocity,  p  is  the  frequency,  and  CQ  is  the 
phase  velocity  at  some  reference  frequency  p  .  Now,  since  the 
frequency  p  is  related  to  the  phase  velocity  C  and  the  wave 
length  A  by: 

_  „  2ttC 
P  -  — 

relation  (a)  can  be  written  as: 

(b)  C  -  Co(l  + 

o 

Here  Aq  is  the  wave  length  at  the  reference  frequency  pQ . 

In  order  to  calculate  tno  value  of  ,  let  us  take  the  total 

a  A 

differentials  of  both  sides  of  expression  (b): 

C  tan  6  AC  tan  6 

dC  =  d[C  +  -2— - An  ~  +  -2 - - - (in  C-£n  A)]  = 

O  71  U 

o 

C  tan  6  C  tan  6  ,  , 

o  dC  o  dA 


Prom  the  last  expression  we  get: 


j 1  tan  6  o>. 

dc<1  — — 


C  tan  6  ,  . 

o  aA 


Hence : 
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For  the  first  harmonic,  this  ratio  has  Its  maximum  absolute 

value  which  is  equal  to  — a~  —  .  Thus,  when  (tan  6)/tt  is  small 

dO 

compared  with  unity,  is  much  smaller  than  C;  for  example, 

dC 

when  tan  6  is  equal  to  0.1,  Ajj-^  is  about  3%  of  the  value  of  C. 


Acknowledgement  s 


92 


The  author  wishes  to  thank  his  adviser,  Professor  H.  Kolsky, 
for  his  kind  assistance. 

Thanks  are  also  due  to  Professor  R.  Clifton  for  helpful 
suggestions,  to  Dr.  J.  Phillips  for  his  help  in  laboratory 
experimentation,  to  Mr.  L.  Daubney  and  Mr.  R.  Stanton  for  their 
careful  preparation  of  the  specimans,  and  to  Mrs.  E.  Fonseca 
for  typing  this  report. 

This  report  was  done  under  a  contract  between  the  Army 
Research  Office  (Durham)  and  Brown  University. 


93 


References 


[1]  B.  Gross,  "Mathematical  Structure  of  the  Theories  of  Visco¬ 
elasticity,"  Hermann,  Paris,  1953. 

[2]  W.  Flflgge,  "Viscoelasticity,"  Blaisdell,  1967. 

[3]  J.D.  Ferry,  "Viscoelastic  Properties  of  Polymers," 

John  Wiley  &  Sons,  New  York,  1961. 

C4]  D.R.  Bland,  "The  Theory  of  Linear  Viscoelasticity," 

Pergamon  Press,  London,  i960. 

[5]  H.  Kolsky,  "Stress  Waves  in  Solids,"  Clarendon  Press, 

Oxford,  1953;  Dover  Reprint,  1964. 

[6]  H.  Kolsky,  "Experimental  Studies  of  the  Mechanical 
Behavior  of  Linear  Viscoelastic  Solids,"  from  "Proc.  of 
the  4th  Symp .  on  Naval  Struc .  Mech . , "  Pergamon  Press, 

Oxford,  New  York,  1966. 

r?]  H.  Kolsky  and  S.S.  Lee,  "The  Propagation  and  Reflection 
of  Stress  Pulses  in  Linear  Viscoelastic  Media,"  Brown 
University  Tech.  Rept .  No.  5,  Contract  Nonr  562(30),  1962. 

[8]  R.A.  Schapery,  "Approximate  Methods  of  Transform  Inversion 
for  Viscoelastic  Stress  Analyses,"  Proc.  4th  U.S.  Nat. 

Cong.  App.  Mech.,  pp.  1075-1085,  1962. 

[9]  K.W.  Hillier,  "A  Method  of  Measuring  some  Dynamic  Elastic 
Constants  and  its  Application  to  the  Study  of  High  Polymers  " 
Proc.  Phys .  Soc . ,  B62,  p.  701,  1949= 

[10]  H.H.  Calvit,  D.  Rader  and  J.  Melville,  "Some  Wave— Propaga¬ 
tion  Experiments  in  Plasteline— clay  Rods,"  Experimental 
Mechanics,  p.  418,  Sept.  1968. 

H.  Kolsky,  "The  Propagation  of  Stress  Pulses  in  Visco¬ 
elastic  Solids,"  Phil.  Mag.,  Vol.  1,  p.  693,  1956. 


[11] 


upper  bar 
lower  bar 
j  unction 


strain  gage 


FIG.  3  DYNAMIC  BRIDGE  CIRCUIT 


/ —  METAL  BLOCK 


/ 


LAPPING  TOOL 


y  sec 
cm 


Fig 


7 


Propagation  of  stress  pulse  from 
aluminum  to  steel. 
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Fig.  11  -  Propagation  of  stress  pulse  from 
polystyrene  to  polyethylene  . 
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Fig.  12a  -  Incident  and  reflected  pulses  in 

polystyrene  bar  attached  to  polyethylene 
bar  (1st  experiment). 
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Fig.  12b  -  Incident  and  reflected  pulses  in 

polystyrene  bar  attached  to  polyethylene 
bar  (2nd  experiment). 


Fig.  14  -  Propagation  of  stress  pulse  from 
aluminum  to  Luc i to  (p.m.m.) 


Pig.  15  -  Incident  and  reflected  pulses  in 

aluminum  bar  attached  to  Lucite  bar 


Fig.  16  -  Recording  of  transmitted  pulse  iru  Lucite 

bur  by  two  strain  gages  mounted  apart. 
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Pig.  17  -  Propagation  of  stress  pulse  from  polystyrene 
to  green  solid  propellant  (g.s.p.) 
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Pig.  18  -  Incident,  and  reflected  pulses  In  polystyrene 
bar  attached  to  g.s.p.  bar 


Fig.  2b  -  Incident  and  reflected  pulses  in  steel 
bar  attached  to  sandstone  bar 


Fig.  2 6  -  Transmitted  pulse  in  sandstone  bar 


Fig.  27  -  Propagation  of  stress  pulse  from 
polystyrene  to  clay 


Fig.  28  -  incident  and  reflected  pulses  in 

polystyrene  har  attached  to  clay  bar 


•i 


i 

i 


Fig.  29  -  Recording  01'  transmitted  pulse  in 

clay  by  two  strain  gages  mounted  8  cm  apart 


Fig.  10  -  Incident  and  reflected  pulses  in 

polystyrene  bar  attached  to  wet  sand  bar 


Fig.  31  -  Incident  and  reflected  pulses  in 

polystyrene  bar  attached  to  dry  sand  bar 
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